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Abstract 


This  document  is  the  Final  Report  deliverable  under  ESTEC  Contract  No.  21303/07/NL/CB 
“GMES  Sentinel- 1  analysis  of  marine  applications  potential  (AMAP)”.  The  document  includes 
ship  detectability  modelling  that  has  been  carried  out  for  Sentinel- 1,  including  development  of  a 
ship  detectability  tool  that  applies  to  Sentinel- 1,  RADARSAT-1,  RADARS AT-2,  and  Envisat 
ASAR  image  data.  Limitations  of  the  tool’s  predictions  are  discussed,  and  examples  from 
RADARSAT-1  and  Envisat  ASAR  are  provided  to  illustrate  the  performance  of  the  tool.  Also 
considered  is  Sentinel- l’s  expected  capability  for  other  maritime  applications,  including  iceberg 
detection,  wind  retrieval,  oil  spill  detection,  sea  ice  surveillance,  and  ship  wake  detection.  In 
general,  the  Interferometric  Wide  Swath  mode  with  its  high  spatial  resolution  and  250  km  swath 
will  be  a  useful  mode  for  maritime  surveillance  in  the  littoral  zone. 


Resume 


Le  present  document  constitue  le  rapport  final  a  presenter  conformement  au  contrat  du  CERS  no 
21303/07/NL/CB  visant  l’analyse,  dans  le  cadre  de  la  mission  Sentinel-1  du  GMES,  du  potentiel 
relatif  a  certaines  applications  maritimes.  II  comprend  une  modelisation  de  la  capacite  de 
detection  de  navires  effectuee  aux  fins  de  la  mission  Sentinel- 1,  dont  la  conception  d’un  outil  de 
prevision  de  la  capacite  de  detection  de  navires  pouvant  utiliser  des  donnees  d’imagerie 
recueillies  par  les  RSO  des  satellites  Sentinel- 1,  RADARSAT-1,  RADARS  AT-2  et  Envisat.  II 
porte  notamment  sur  les  limites  des  previsions  de  V  outil  et  comprend  des  exemples  obtenus  au 
moyen  des  RSO  de  RADARSAT-1  et  d’Envisat,  afin  d’illustrer  le  rendement  de  V outil.  On  y 
traite  aussi  de  la  capacite  prevue  de  Sentinel- 1,  en  ce  qui  concerne  d’autres  applications 
maritimes,  dont  la  detection  d’ icebergs,  V extraction  des  vents,  la  detection  des  deversements  de 
petrole,  la  surveillance  des  glaces  de  mer  et  la  detection  du  sillage  de  navires.  En  general,  le  mode 
d’interferometrie  a  fauchee  large,  dont  la  resolution  spatiale  est  elevee  et  la  fauchee  mesure  250 
km,  se  prete  a  la  surveillance  maritime  dans  les  zones  littorales. 
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Executive  summary 


GMES  Sentinel-1  Analysis  of  Marine  Applications  Potential 
(AMAP) 

Paris  W.  Vachon;  John  Wolfe;  DRDC  Ottawa  ECR  2008-218;  Defence  R&D 
Canada  -  Ottawa;  October  2008. 

Introduction:  The  ESA  Sentinels  are  operational  satellites  that  will  support  the  Global 
Monitoring  for  Environment  and  Security  (GMES)  program.  The  Sentinel- 1  mission  will  consist 
of  a  constellation  of  two  satellites  that  will  carry  C-band  synthetic  aperture  radar  (SAR).  These 
radars  will  support  a  variety  of  applications  including  maritime  surveillance.  Whereas  Sentinel- 1 
will  have  several  modes  of  operation,  the  focus  of  this  work  has  been  on  the  Interferometric  Wide 
(IW)  mode  with  its  250  km  swath,  5x20  metre  spatial  resolution  (ground-range  by  azimuth,  worst 
case),  a  single  look,  and  selectable  dual  polarization  (i.e.,  HH+HV  or  VV+VH).  IW  is  expected 
to  be  the  default  mode  of  operation  in  the  littoral  zone.  Performance  of  the  IW  mode  for  maritime 
surveillance  applications,  including  ship  detection  iceberg  detection,  wind  retrieval,  oil  spill 
detection,  sea  ice  detection,  and  ship  wake  detection  has  been  compared  to  current  operational 
practice  with  RADARSAT-1  and  Envisat  ASAR. 

Results:  The  focus  of  this  work  is  on  ship  detectability  modelling,  which  has  been  carried  out  for 
Sentinel- 1,  including  development  of  a  ship  detectability  performance  tool  that  applies  to 
Sentinel- 1,  RADARSAT-1,  RADARSAT-2,  and  Envisat  ASAR  image  data.  The  tool  assumes 
straightforward,  pixel-based  detection  of  ship  targets  in  background  ocean  clutter  and  includes 
separate  elements  describing  the  ship  radar  cross  section,  the  ocean  clutter,  and  the  radar  system 
noise.  Limitations  of  the  tool,  including  the  impacts  of  target  motion,  are  discussed,  and  are 
included  in  the  tool  through  a  user-selectable  margin  parameter.  Examples  from  RADARSAT-1 
and  Envisat  ASAR  are  provided  to  demonstrate  the  performance  of  the  tool. 

Ship  detectability  modelling  has  shown  that,  due  to  its  spatial  resolution  and  range  of  incidence 
angles,  IW  mode  is  very  well-suited  to  ship  detection,  almost  regardless  of  the  available 
polarization.  HH  polarization  is  the  best  of  the  co-polarization  choices,  with  cross-polarization 
improving  performance  at  smaller  incidence  angles.  Operational  use  of  two  polarizations  would 
improve  the  detection  performance  or  could  be  used  to  reject  false  alarms. 

Other  maritime  applications  were  also  evaluated  based  upon  backscatter  models  and  imagery 
from  RADARSAT-1  or  Envisat  ASAR.  It  is  apparent  that  IW  will  be  well  suited  to:  iceberg 
detection  and  ship/iceberg  discrimination;  wind  retrieval,  though  smaller  incidence  angles  would 
be  better;  oil  spill  detection,  though  smaller  incidence  angles  would  be  better;  sea  ice 
surveillance,  though  a  larger  swath  would  be  desirable;  and  ship  wake  detection,  though  smaller 
incidence  angles  would  be  better. 

As  a  general  statement,  the  one-size-fits-all  approach  that  will  be  adopted  with  Sentinel- 1  IW 
operations  will  suit  most  maritime  applications  in  the  littoral  zone  rather  well. 

Significance:  This  document  is  the  Final  Report  deliverable  under  ESTEC  Contract  No. 
21303/07/NL/CB  “GMES  Sentinel- 1  analysis  of  marine  applications  potential  (AMAP)”. 
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GMES  Sentinel-1  Analysis  of  Marine  Applications  Potential 
(AMAP) 

Paris  W.  Vachon;  John  Wolfe;  DRDC  Ottawa  ECR  2008-218;  R  &  D  pour  la 
defense  Canada  -  Ottawa;  Octobre  2008. 

Introduction  :  Les  satellites  Sentinel  de  l’ASE  sont  des  engins  operationnels  qui  appuieront  le 
programme  GMES  (Global  Monitoring  for  Environment  and  Security  -  surveillance  mondiale 
pour  l’environnement  et  la  securite).  La  mission  Sentinel-1  est  fondee  sur  une  constellation  de 
deux  satellites  dotes  de  radars  a  synthese  d’ouverture  (RSO)  en  bande  C,  lesquels  appuieront 
diverses  applications,  dont  la  surveillance  maritime.  Bien  que  la  mission  Sentinel- 1  implique 
plusieurs  modes  de  fonctionnement,  elle  est  axee  sur  le  mode  d’interferometrie  a  fauchee  large, 
d’une  fauchee  de  250  km,  d’une  resolution  spatiale  de  5  m  sur  20  m  (distance  au  sol  d’apres 
l’azimut,  dans  le  pire  des  cas),  a  visee  unique  et  a  polarisation  double  (soit  HH+HV  ou  VV+VH). 
Ce  mode  d’interferometrie  devrait  etre  le  mode  de  fonctionnement  par  defaut  dans  la  zone 
littorale.  Son  rendement  a  des  fins  de  surveillance  maritime,  dont  la  detection  d’icebergs, 
1’ extraction  des  vents,  la  detection  des  deversements  de  petrole,  la  detection  des  glaces  de  mer  et 
la  detection  du  sillage  de  navires,  a  ete  compare  a  celui  des  modes  de  fonctionnement  courants 
avec  les  RSO  de  RADARSAT-1  et  d’Envisat. 

Resultats  :  Les  presents  travaux  sont  axes  sur  la  modelisation  de  la  capacite  de  detection  de 
navires  effectuee  dans  le  cadre  de  la  mission  Sentinel- 1,  dont  la  conception  d’un  outil  de 
prevision  de  la  capacite  de  detection  de  navires  pouvant  utiliser  des  donnees  d’imagerie 
recueillies  par  les  RSO  des  satellites  Sentinel- 1,  RADARSAT-1,  RADARS AT-2  et  Envisat. 
L’outil  se  base  sur  la  detection  simple  par  pixels  de  navires  cibles  au  sein  d’un  fouillis  d’echos 
dans  1’ ocean  et  comprend  des  elements  distincts  decrivant  la  section  efficace  en  radar  des  navires, 
le  fouillis  d’echos  dans  T ocean  et  le  bruit  du  radar.  Les  limites  de  l’outil,  dont  l’incidence  du 
deplacement  des  cibles,  sont  abordees  dans  le  present  document  et  prises  en  consideration  par 
l’outil,  par  le  biais  d’un  parametre  de  marge  que  l’utilisateur  peut  choisir.  Des  exemples  obtenus 
au  moyen  des  RSO  de  RADARSAT-1  et  d’Envisat  sont  presentes,  afin  de  demontrer  le 
rendement  de  l’outil. 

La  modelisation  de  la  capacite  de  detection  de  navires  a  montre  qu’en  raison  de  sa  resolution 
spatiale  et  de  son  eventail  d’ angles  d’ incidence,  le  mode  d’interferometrie  a  fauchee  large  se  prete 
bien  a  la  detection  de  navires,  peu  importe  la  polarisation  disponible,  ou  presque.  La  polarisation 
double  HH  est  la  meilleure,  la  polarisation  croisee  ameliorant  la  performance  a  des  angles 
d’incidence  inferieurs.  L’utilisation  operationnelle  de  deux  polarisations  pourrait  ameliorer  la 
capacite  de  detection  ou  servir  a  eviter  de  fausses  alarmes. 

D’autres  applications  maritimes  ont  ete  evaluees  d’apres  des  modeles  de  retrodiffusion  et  des 
images  produites  par  les  RSO  de  RADARSAT-1  ou  d’Envisat.  L’interferometrie  a  fauchee  large 
se  pretera  manifestement  bien  a  la  detection  d’icebergs,  a  la  distinction  des  icebergs  et  des 
navires,  a  l’extraction  des  vents  (des  angles  d’incidence  inferieurs  seraient  preferables),  a  la 
detection  des  deversements  de  petrole  (des  angles  d’incidence  inferieurs  seraient  preferables),  a  la 
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detection  des  glaces  de  mer  (une  fauchee  plus  large  serait  preferable)  et  a  la  detection  du  sillage 
de  navires  (des  angles  d’incidence  inferieurs  seraient  preferables). 

En  general,  l’approche  uniformisee  qui  sera  utilisee  pendant  les  operations  d’interferometrie  a 
fauchee  large  de  la  mission  Sentinel- 1  se  pretera  assez  bien  a  la  plupart  des  applications 
maritimes  dans  la  zone  littorale. 

Importance  :  Le  present  document  constitue  le  rapport  final  a  presenter  conformement  au  contrat 
du  CERS  no  21303/07/NL/CB  visant  l’analyse,  dans  le  cadre  de  la  mission  Sentinel-1  du  GMES, 
du  potentiel  relatif  a  certaines  applications  maritimes. 
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1  Introduction 


The  ESA  Sentinels  are  operational  satellites  that  will  support  the  Global  Monitoring  for 
Environment  and  Security  (GMES)  program  [1].  The  Sentinel- 1  mission  will  consist  of  a 
constellation  of  two  satellites  that  will  carry  C-band  synthetic  aperture  radar  (SAR).  These  radars 
will  support  a  variety  of  applications  including  maritime  surveillance. 

A  key  question  is:  “What  is  the  capability  of  Sentinel- 1  for  maritime  surveillance  and  ship 
detection  in  particular?”  In  this  project  “GMES  Sentinel- 1  analysis  of  marine  applications 
potential  (AMAP)”,  this  question  has  been  addressed  by  carrying  out  the  following  tasks: 

1.  Create  a  model  for  ship  detection  performance  prediction  based  upon  Vachon  et  al.  (1997) 
[30]  and  extend  the  model  based  upon  recent  ship  detection  work  by  Vachon  et  al.  [32]  [34] 
[36]  that  also  addresses  ship  detectability  by  the  Envisat  ASAR  AP  mode  of  operation; 

2.  Implement  the  model  in  Matlab®  such  that  the  performance  of  all  Sentinel- 1  modes  are 
addressed  and  may  be  compared  with  the  performance  of  RADARSAT-1,  RADARS AT-2, 
and  Envisat  ASAR; 

3.  Develop  demonstration  products  using  RADARSAT-1  and  Envisat  ASAR  data  that  are 
currently  available  at  DRDC  Ottawa;  and 

4.  Address  the  capability  of  the  Interferometric  Wide  (IW)  mode  and  the  other  Sentinel- 1  modes 
for  iceberg  detection,  wind  retrieval,  oil  spill  detection,  sea  ice  detection,  and  ship  wake 
detection  based  upon  the  expected  clutter-to-noise  ratio,  wind  speed  sensitivity,  oil  and  ice 
backscatter,  and  available  observational  evidence. 

This  document  is  the  Final  Report  project  deliverable;  it  addresses  items  1  through  4  of  the  above 
tasks. 

1.1  Maritime  surveillance  requirements 

There  is  a  vast  literature  that  addresses  the  role  of  SAR  in  the  ocean  environment,  e.g.,  [7],  [12]. 

For  maritime  surveillance  applications,  wide  swath  coverage  is  generally  required  to  maximize 
the  rate  of  area  coverage.  This  requirement  is  fundamental. 

Ship  detection  relies  on  the  detection  of  bright  point-like  targets  against  an  ocean  clutter 
background.  The  ocean  clutter  level  is  influenced  predominantly  by  the  local  wind  speed,  which 
generates  centimetre-scale  Bragg  scattering  waves  or  “roughness”  on  the  ocean  surface.  The 
higher  the  wind,  the  rougher  the  surface,  and  the  higher  the  ocean  clutter  level.  In  fact, 
radiometrically  calibrated  radar  can  be  used  to  retrieve  high  resolution  wind  fields.  However,  the 
clutter  is  also  subject  to  local  modulations  from  various  physical  processes,  such  as  local  wind 
speed  variations,  wakes,  and  atmospheric  turbulence.  Oceanographic  modulations  due  to 
(oceanic)  internal  waves,  current  shear,  and  the  presence  of  surfactants  that  damp  the  surface 
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roughness,  also  affect  the  clutter.  Ship  detection  is  most  effective  when  the  ocean  clutter  level  is 
low.  The  clutter  decreases  with  increasing  incidence  angle. 

Iceberg  detection  is  closely  related  to  ship  detection.  Iceberg  signatures  tend  to  have  different 
shapes  and  wake  characteristics  than  those  of  ships,  and  different  signatures  as  a  function  of  radar 
polarization. 

Oil  spill  detection  relies  upon  the  suppression  of  clutter  by  the  presence  of  oil.  This  means  that 
the  background  clutter  must  be  sufficiently  above  the  instrument  noise  floor  to  allow  clutter 
suppression  to  be  observed.  As  the  wind  speed  increases,  the  oil  layer  breaks  up  and  the  oil 
mixes  downwards  into  the  water  column.  Detection  of  ships  in  the  vicinity  of  detected  oil  spills 
provides  information  that  may  be  helpful  to  associate  a  polluter  with  a  particular  oil  spill. 

Sea  ice  surveillance  requires  measurement  of  the  ice  edge  and  classification  of  the  ice  type: 
smooth  new,  first  year,  and  multi-year  ice.  The  ice  edge  is  more  easily  detected  when  the  sea/ice 
contrast  is  large,  that  is,  when  the  ocean  clutter  level  is  low,  as  occurs  for  larger  incidence  angles. 
Discrimination  between  ice  types  is  also  incidence  angle  dependent,  with  discrimination  between 
ice  types  increasing  with  increasing  incidence  angle. 

Observation  of  ship  wakes  is  also  of  interest.  The  presence  of  a  ship  wake  modulates  the 
background  clutter.  Measurement  of  ship  wake  properties  can  allow  estimation  of  the  ship  speed 
and  course.  Ship  wakes  are  easier  to  detect  if  the  clutter  level  is  high,  although  the  wake 
signature  may  be  obliterated  if  the  wind  speed  is  too  large. 

There  are  several  competing  factors  when  trying  to  meet  a  maritime  surveillance  requirement 
with  SAR.  The  key  radar  parameters  are: 

•  Swath  width; 

•  Incidence  angle; 

•  Spatial  resolution; 

•  Radiometric  resolution;  and 

•  Noise  floor. 

In  this  project,  straightforward  C-band  radar  backscatter  models  along  with  nominal  radar 
parameters  are  used  to  evaluate  the  expected  performance  of  Sentinel- 1  for  ship  detection,  wind 
retrieval,  oil  spill  detection,  sea  ice  surveillance,  and  wake  detection. 

1.2  Outline  of  document 

Section  2  of  this  document  includes  a  description  of  the  Sentinel- 1  modes  and  the  nominal  sensor 
characteristics  that  have  been  used.  Also  included  are  the  nominal  parameters  of  RADARS  AT- 1, 
RADARS AT-2,  and  Envisat  ASAR,  for  comparison  purposes.  Section  3  considers  ship 
detectability  modelling,  which  leads  to  a  minimum  detectable  ship  length  that  permits 
intercomparison  of  various  radars  and  modes.  A  parametric  study  shows  the  model  sensitivity 
and  several  possible  system  design  trade-offs.  Section  4  focuses  on  iceberg  detection,  wind 
retrieval,  oil  spill  detection,  sea  ice  detection,  and  ship  wake  detection.  Conclusions  are  offered 
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in  Section  5.  The  Annexes  contain  a  description  of  the  ship  detectability  model  as  implemented 
in  Matlab®  and  some  model  runs  along  with  RADARSAT-1  and  Envisat  ASAR  imagery  of 
known  ships. 
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2  Sentinel-1  description 


Sentinel- 1  will  offer  four  modes  of  operation: 

1.  Strip  Map  (S)  Mode  offers  six  different  beams  with  an  80  km  swath,  5x5  metre  spatial 
resolution  (ground-range  by  azimuth,  worst  case),  a  single  look,  and  selectable  dual 
polarization  (i.e.,  HH+HV  or  VV+VH); 

2.  Interferometric  Wide  (IW)  Swath  Mode  is  composed  of  three  different  beams  providing  a  250 
km  swath,  5x20  metre  spatial  resolution  (ground-range  by  azimuth,  worst  case),  a  single  look, 
and  selectable  dual  polarization  (i.e.,  HH+HV  or  VV+VH); 

3.  Extra-Wide  (EW)  Swath  Mode  is  composed  of  five  different  beams  providing  a  400  km 
swath,  20x40  metre  spatial  resolution  (ground-range  by  azimuth,  worst  case),  a  single  look, 
and  selectable  dual  polarization  (i.e.,  HH+HV  or  VV+VH);  and 

4.  Wave  (W)  Mode  with  5x20  metre  spatial  resolution  (ground-range  by  azimuth,  worst  case),  a 
single  look,  selectable  single  polarization,  and  sampled  images  of  20x20  km  every  100  km 
along  the  orbit  track. 

The  Sentinel- 1  mission  philosophy  is,  to  the  extent  possible,  to  operate  continuously  in  a  single 
mode  over  land  and  the  coastal  zone.  The  preferred  mode  is  IW,  with  Wave  Mode  being  used  for 
open  ocean  areas.  The  focus  of  this  work  is  on  IW  performance,  with  comments  on  the  S  and  EW 
modes  where  applicable. 

Sentinel- 1  will  be  operated  in  a  circular  orbit  with  a  reference  altitude  of  711  km  relative  to  the 
WGS84  ellipsoid.  However,  there  will  be  a  slight  orbit  eccentricity  with  respect  to  the  Earth 
centre;  the  satellite  altitude  will  vary  along  the  orbit  by  about  28  km.  In  order  to  minimize  the 
number  of  pulse  repetition  frequency  (PRF)  changes,  a  roll  steering  law  has  been  defined.  The 
roll  steering  will  reduce  the  slant-range  variation  along  the  orbit,  leading  to  a  changing  incidence 
angle  relative  to  a  fixed  point  in  the  swath.  In  this  work,  the  Sentinel- 1  parameters  at  the  extreme 
altitude  points  along  the  orbit,  namely  the  low-orbit  (697  km)  and  high-orbit  (725  km)  cases,  have 
been  modelled.  The  actual  performance  achieved  at  some  other  point  along  the  orbit  would  lie 
between  these  two  extremes. 

Nominal  mid-swath  values  for  the  relevant  Sentinel- 1  parameters  are  summarized  in  Table  1  in 
which  0Min  is  the  minimum  incidence  angle  within  the  beam;  dMax  is  the  maximum  incidence  angle 
within  the  beam;  pA  is  the  azimuth  resolution;  pSR  is  the  slant-range  resolution;  ENL  is  the 
equivalent  number  of  independent  looks;  and  NESZ  is  the  noise-equivalent  sigma-zero  (worst 
case  within  the  beam).  For  beam-oriented  calculations,  the  actual  noise  floor  and  ground  range 
resolution  (both  functions  of  the  local  incidence  angle)  are  used.  Detailed  Sentinel- 1  information 
has  been  compiled  from  available  literature,  including  [1],  [23],  and  [24],  along  with  input  from 
the  ESA  Technical  Authority  for  this  project. 
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Table  1:  Nominal  Sentinel- 1  parameters. 


Low  orbit  (697  km) 

High  orbit  (725  km) 

Beam 

® Min 

@ Max 

Pa 

PSR 

ENL 

NESZ 

P  Min 

6  Max 

Pa 

PSR 

ENL 

NESZ 

[deg] 

[deg] 

[m] 

[m] 

[dB] 

[deg] 

[deg] 

[m] 

[m] 

[dB] 

SI 

20.0 

26.3 

4.9 

1.5 

1 

-24.0 

18.3 

24.5 

4.9 

1.5 

1 

-24.4 

S2 

23.5 

29.5 

4.9 

1.8 

1 

-23.4 

21.8 

27.7 

4.9 

1.8 

1 

-23.6 

S3 

29.3 

34.8 

4.9 

2.3 

1 

-25.0 

27.6 

33.0 

4.9 

2.3 

1 

-24.7 

S4 

34.7 

39.7 

4.9 

2.7 

1 

-23.6 

33.0 

38.0 

4.9 

2.7 

1 

-23.7 

S5 

39.6 

44.1 

4.9 

3.0 

1 

-24.2 

37.9 

42.4 

4.9 

3.0 

1 

-23.4 

S6 

42.5 

46.7 

4.9 

3.2 

1 

-24.4 

40.8 

45.0 

4.9 

3.2 

1 

-23.7 

IW1 

30.9 

36.6 

19.5 

2.4 

1 

-24.1 

29.2 

34.9 

19.5 

2.4 

1 

-23.8 

IW2 

36.5 

41.8 

19.5 

2.8 

1 

-23.6 

34.8 

40.1 

19.5 

2.8 

1 

-23.2 

IW3 

41.7 

46.0 

19.5 

3.2 

1 

-23.6 

40.0 

44.2 

19.5 

3.2 

1 

-23.1 

EW1 

20.8 

29.3 

40.0 

6.5 

1 

-25.8 

19.0 

27.5 

40.0 

6.5 

1 

-25.8 

EW2 

29.2 

35.1 

40.0 

9.2 

1 

-27.2 

27.4 

33.3 

40.0 

9.2 

1 

-27.2 

EW3 

35.1 

39.8 

40.0 

11.0 

1 

-27.2 

33.3 

38.0 

40.0 

11.0 

1 

-27.2 

EW4 

39.7 

44.2 

40.0 

12.3 

1 

-25.6 

38.0 

42.5 

40.0 

12.3 

1 

-25.6 

EW5 

44.3 

47.0 

40.0 

13.5 

1 

-32.3 

42.5 

45.2 

40.0 

13.5 

1 

-32.2 

2.1  Other  sensors  of  interest 

Sentinel-1  performance  will  be  considered  in  the  context  of  RADARSAT-1,  RADARSAT-2,  and 
Envisat  ASAR,  with  the  relevant  parameters  for  these  sensors  listed  in  Table  2  through  Table  4. 
The  SAR  system  information  was  derived  from  [5],  [22],  and  [26]. 

RADARSAT-1  is  an  HH  polarization  instrument.  As  for  Sentinel- 1,  nominal  parameters  (Table 
2)  are  provided  and  the  indicated  noise  floors  are  the  worst  case  within  the  swath,  unless  written 
as  “-25”  in  which  case  they  are  a  nominal  value.  RADARSAT-1  has  two  ScanSAR  modes 
(Table  3):  ScanSAR  Narrow  (SCN)  is  composed  of  2  (SCNA)  or  3  (SCNB)  beams  with  a 
nominal  300  km  swath;  ScanSAR  Wide  (SCW)  is  composed  of  3  (SCWA)  or  4  (SCWB)  beams 
with  a  nominal  500  km  swath. 

RADARSAT-2  data  are  modelled  with  the  RADARSAT-1  parameters  for  common  (i.e.,  heritage) 
modes,  but  with  the  following  polarization  choices  available  for  RADARSAT-2:  selective  single 
polarization  (HH  or  VV)  or  dual  polarization  (HH+HV  or  VV+VH).  The  RADARSAT-2  noise 
floor  is  taken  as  5  dB  smaller  than  that  of  RADARSAT-1.  The  higher  resolution  and  polarimetric 
modes  offered  by  RADARSAT-2  have  not  been  considered  in  this  work. 

Envisat  ASAR  offers  selective  single  polarization  (IMP  product),  alternating  polarization  (APP 
product)  that  trades  radiometric  resolution  for  polarization  in  a  burst-mode  of  operation,  and  a 
wide  swath  mode  (WSM)  composed  of  5  beams  with  a  nominal  swath  width  of  450  km  (Table  4). 
For  APP  mode,  the  ground  range  resolution  is  traded  for  ENL  for  larger  incidence  angle  beams. 
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Table  2:  Nominal  RADARSAT-1  single  beam  parameters. 


Beam 

0 Min 

[deg] 

0Max 

[deg] 

Pa 

[m] 

PSR 

[m] 

ENL 

NESZ 

[dB] 

SI 

19.5 

26.7 

27.0 

9.3 

3.1 

-25.7 

S2 

24.1 

30.9 

27.0 

9.3 

3.1 

-24.4 

S3 

31.0 

37.0 

27.0 

14.1 

3.1 

-26.5 

S4 

33.6 

39.4 

27.0 

14.1 

3.1 

-26.4 

S5 

36.4 

41.9 

27.0 

14.1 

3.1 

-25.8 

S6 

41.7 

46.5 

27.0 

14.1 

3.1 

-25.7 

S7 

44.7 

49.2 

27.0 

14.1 

3.1 

-25.8 

FI 

36.8 

39.9 

8.4 

5.2 

1 

-25.0 

F2 

39.3 

42.1 

8.4 

5.2 

1 

-23.9 

F3 

41.5 

44.1 

8.4 

5.2 

1 

-25.4 

F4 

43.5 

45.8 

8.4 

5.2 

1 

-25.8 

F5 

45.4 

47.6 

8.4 

5.2 

1 

-25.0 

W1 

19.3 

30.2 

27.0 

14.1 

3.1 

-24.3 

W2 

30.2 

38.9 

27.0 

14.1 

3.1 

-23.9 

W3 

48.9 

45.1 

27.0 

14.1 

3.1 

-23.6 

EH1 

49.0 

53.4 

27.0 

14.1 

3.1 

-25 

EH2 

50.1 

53.5 

27.0 

14.1 

3.1 

-25 

EH3 

51.2 

54.6 

27.0 

14.1 

3.1 

-25 

EH4 

54.5 

57.2 

27.0 

14.1 

3.1 

-25 

EH5 

55.6 

58.2 

27.0 

14.1 

3.1 

-25 

EH6 

56.9 

59.4 

27.0 

14.1 

3.1 

-25 

ELI 

10.4 

21.9 

27.0 

9.3 

3.1 

-24.5 
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Table  3:  Nominal  RADARSAT-1  ScanSAR parameters. 


Mode 

ft  Min 

Tdegl 

ft  Max 

fdegl 

Pa 

[ml 

P  SR 

rml 

ENL 

NESZ 

[dBl 

SCNA 

W1 

19.3 

30.2 

47.8 

27.1 

3.5 

-24.3 

W2 

30.2 

38.9 

53.8 

27.1 

3.5 

-23.9 

SCNB 

W2 

30.2 

36.4 

71.1 

27.1 

3.5 

-23.9 

S5 

36.4 

41.7 

71.9 

27.1 

3.5 

-25.8 

S6 

41.7 

46.9 

78.8 

27.1 

3.5 

-25.7 

SCWA 

W1 

19.3 

30.2 

93.1 

53.8 

7.0 

-24.3 

W2 

30.2 

38.9 

104.7 

53.8 

7.0 

-23.9 

W3 

38.9 

45.1 

117.3 

53.8 

7.0 

-23.6 

S7 

45.1 

49.2 

117.5 

53.8 

7.0 

-25.8 

SCWB 

W1 

19.3 

30.2 

93.1 

53.8 

7.0 

-24.3 

W2 

30.2 

36.4 

104.7 

53.8 

7.0 

-23.9 

S5 

36.4 

41.7 

106.0 

53.8 

7.0 

-25.8 

S6 

41.7 

46.9 

117.6 

53.8 

7.0 

-25.7 

Table  4:  Nominal  Envisat  ASAR  parameters. 


Beam 

ft  Min 

[deg] 

ft  Max 

[deg] 

Pa 

[m] 

PSR 

[m] 

ENL 

NESZ 

rdBl 

IMP-IS1 

15.0 

22.9 

22.1 

9.6 

3.95 

-19.6 

IMP-IS2 

19.2 

26.7 

22.1 

9.4 

3.95 

-19.6 

IMP-IS3 

26.0 

31.4 

22.1 

11.8 

3.95 

-20.6 

IMP-IS4 

31.0 

36.3 

22.1 

14.2 

3.95 

-19.0 

IMP-IS5 

35.8 

39.4 

22.1 

15.8 

3.95 

-19.0 

IMP-IS6 

39.1 

42.8 

22.1 

17.2 

3.95 

-21.2 

IMP-IS7 

42.5 

45.2 

22.1 

18.5 

3.95 

-20.7 

APP-IS1 

15.0 

22.9 

27.7 

9.7 

1.76 

-19.6 

APP-IS2 

19.2 

26.7 

27.7 

9.6 

1.73 

-19.6 

APP-IS3 

26.0 

31.4 

27.7 

12.1 

2.25 

-20.6 

APP-IS4 

31.0 

36.3 

27.7 

14.3 

2.66 

-19.0 

APP-IS5 

35.8 

39.4 

27.7 

16.0 

3.30 

-19.0 

APP-IS6 

39.1 

42.8 

27.7 

17.2 

3.78 

-21.2 

APP-IS7 

42.5 

45.2 

27.7 

18.4 

3.73 

-20.7 

SSI 

17.0 

26.0 

107.5 

44.9 

13.19 

-19.9 

SS2 

26.0 

31.5 

107.5 

54.4 

13.21 

-23.4 

SS3 

31.0 

36.5 

107.5 

63.6 

13.84 

-21.8 

SS4 

36.0 

39.0 

107.5 

73.6 

13.77 

-23.1 

SS5 

38.5 

42.6 

107.5 

78.7 

13.38 

-25.4 
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3  Ship  detectability 


3.1  Problem  description 

In  SAR  imagery,  ships  are  generally  detected  by  using  pixel-based  algorithms  that  search  for 
bright  pixels  against  the  ocean  clutter  background.  Ship  detection  is  most  effective  when  the 
ocean  clutter  level  is  low  and  the  spatial  resolution  of  the  sensor  is  high.  Modelling  ship 
detectability  is  a  statistical  problem  that  considers  separate  models  for  the  system  noise,  the  ocean 
clutter,  and  the  ship  radar  cross  section  (RCS)  and  its  variability  (i.e.,  target  fading)  [30].  The 
detection  problem  is  illustrated  in  Figure  1.  The  two  distributions  represent  the  clutter  plus  noise 
(on  the  left)  and  the  ship  RCS  variability  (on  the  right).  The  detection  performance  is  generally 
parameterized  by  the  probability  of  false  alarm  (Pfa)  and  the  probability  of  detection  ( PD ).  The 
area  under  the  clutter  plus  noise  distribution  up  to  the  detection  threshold  line  is  (1  -  Pfa)-  The 
area  under  the  ship  RCS  variability  distribution  down  to  the  detection  threshold  line  is  PD.  A  key 
requirement  in  a  practical  ship  detection  system  is  to  choose  the  detection  threshold  such  that  PFA 
is  extremely  small  (e.g.,  <  2.5(10)'9)  and  PD  is  large  (e.g.,  >  0.8).  It  is  generally  helpful  to  express 
the  ship  detection  performance  in  terms  of  a  ship  length  rather  than  a  radar  cross  section, 
requiring  a  model  for  the  relationship  between  ship  RCS  and  ship  length. 


Figure  1:  Illustration  of  the  ship  detection  problem  in  which  the  clutter  PDF  (blue)  and  the  ship 
PDF  (red)  overlap,  with  the  solid  black  line  representing  the  threshold  that  sets  PFa  and  PD. 
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3.2  Ship  detection  model 

In  this  section,  a  model  for  ship  detectability  is  developed. 

3.2.1  System  noise 

The  system  noise  arises  in  the  receiver  and  is  modelled  as  additive  Gaussian  noise.  However,  it  is 
generally  expressed  as  a  noise-equivalent  normalized  RCS,  or  NESZ  (noise-equivalent  sigma- 
zero),  which  depends  upon  the  range  to  the  scene,  the  incidence  angle,  and  the  elevation  antenna 
pattern,  among  other  things.  NESZ  is  a  useful  concept  since  it  can  be  compared  directly  with  the 
observed  backscatter  level  to  establish  the  clutter-to-noise  ratio  (CNR).  An  example  of  an  NESZ 
plot  is  shown  in  Figure  2.  The  NESZ  is  variable  with  range  (lumpy  in  this  example)  due  mainly 
to  the  shape  of  the  elevation  antenna  pattern. 

The  system  noise  would  be  inverse-exponentially  distributed  for  a  single-look  SAR  system. 
Then,  for  a  multi-look  system,  which  may  be  characterized  by  its  equivalent  number  of 
statistically  independent  looks  (ENL),  the  noise  becomes  chi-squared  distributed  with  2*ENL 
degrees-of- freedom. 


Figure  2:  Plots  of  ocean  clutter  level  as  a  function  of polarization ,  wind  speed ,  and  wind 
direction,  along  with  a  plot  of  the  noise  floor  for  a  RADARSAT-1  SCW  product. 
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3.2.2  Ocean  clutter 


The  ocean  clutter  level  is  influenced  predominantly  by  the  local  wind  speed,  which  generates 
centimetre-scale  Bragg  scattering  waves  or  “roughness”  on  the  ocean  surface.  The  higher  the 
wind,  the  rougher  the  surface,  and  the  higher  the  ocean  clutter  level.  However,  the  clutter  is  also 
subject  to  local  modulations  from  various  physical  processes,  such  as  local  wind  speed  variations, 
wakes,  and  turbulence.  Oceanographic  modulations  due  to  (oceanic)  internal  waves,  current 
shear,  and  the  presence  of  surfactants,  which  damp  the  roughness,  also  affect  the  clutter. 

The  normalized  RCS  (cr°)  or  backscatter  is  a  measure  of  the  surface  roughness.  The  ocean  surface 
roughness  is  governed  primarily  by  the  local  wind  conditions.  The  normalized  RCS  can  be  used 
to  estimate  the  wind  speed  at  the  ocean  surface  to  high  spatial  resolution  (1  km  or  better,  as 
compared  to  resolutions  on  the  order  of  10  km,  for  a  spaceborne  scatterometer).  By  using 
principles  from  wind  scatterometry,  a  semi-empirical  model  function  may  be  used  to  estimate  the 
surface  wind  speed  from  the  observed  normalized  RCS,  the  local  incidence  angle,  and  the  relative 
wind  direction.  The  normalized  RCS  may  be  estimated  directly  from  the  SAR  image;  the  local 
incidence  angle  is  derived  from  the  acquisition  geometry.  It  remains  a  challenge  to  estimate  the 
wind  direction  from  SAR  image  data  alone,  but  there  are  often  useful  indicators  in  the  imagery 
such  as  boundary  layer  rolls,  wind  shadows,  and  wakes,  or  external  sources  of  wind  direction 
such  as  numerical  weather  prediction  models,  buoy  observations,  and  scatterometers  may  be  used. 

In  this  work,  the  following  models  for  the  C-band  normalized  RCS  (er°)  of  the  ocean  surface  are 
considered.  For  VV  polarization  CMOD-IFREMER2,  a  C-band  wind  scatterometry  model  that 
was  developed  at  IFREMER  (France)  for  processing  of  archived  ERS  scatterometer  data  is  used; 
for  HH  polarization  a  hybrid  model  composed  of  CMOD-IFREMER2  and  a  co-polarization  ratio 
that  is  based  upon  Kirchhoff  scattering  is  used.  These  models  have  been  extensively  discussed 
and  validated  in  the  RADARSAT-1  context,  e.g.,  [31].  Representative  clutter  transects  based 
upon  this  model  are  also  shown  in  Figure  2. 

Expressed  in  decibels,  the  clutter-to-noise  ratio  (CNR)  for  a  particular  wind  speed,  direction,  and 
incidence  angle  is  the  difference  between  a  particular  normalized  RCS  curve  and  the  noise  floor. 
The  CNR  approaches  zero  decibels  (i.e.,  equal  clutter  and  noise  contributions)  for  increasing 
incidence  angles.  Negative  clutter-to-noise  ratios  indicate  noise-limited  imaging  (i.e.,  the  SAR 
image  is  dominated  by  system  noise). 

The  CNR  increases  for  increasing  wind  speed.  For  the  same  wind  speed  and  direction,  VV 
polarization  will  always  have  a  larger  CNR  than  HH  polarization.  For  the  same  wind  speed  and 
polarization,  an  upwind  or  a  downwind  look  direction  (i.e.,  wind  blowing  towards  or  away  from 
the  radar)  will  always  have  a  larger  CNR  than  an  across-wind  look  direction  (i.e.,  wind  blowing 
perpendicular  to  the  radar  look  direction).  The  wind  direction  sensitivity  is  caused  by  the 
directional  distribution  of  the  Bragg-scale  waves. 

Cross  polarization  requires  a  further  modification  to  the  model.  It  is  well-known  that  cross-pol 
clutter  is  extremely  small  and  is  generally  expected  to  be  near  or  below  the  noise  floor  of  the 
instrument.  An  airborne  SAR  example  is  presented  in  [33].  For  many  spaceborne  systems,  it 
may  be  assumed  that  the  cross-pol  clutter  level  is  zero  such  that  ship  detection  is  against  the 
instrument  noise  floor  rather  than  against  the  background  clutter.  Recent  work  with  Envisat 
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ASAR  AP  mode  data  shows  that  there  could  be  some  cross-pol  clutter  signature  at  smaller 
incidence  angles  [36]. 

Based  upon  fitting  to  the  data  presented  in  [33],  a  very  simple  formulation  for  the  cross-pol  clutter 
has  been  developed  as  a  correction  to  the  VV  polarization  clutter: 

<j°vh  =  cr°HV  =  <r°^ (£/,#,<£>) +  0.25  U  +  0.22  0  -  28.1  [dB]  (1) 

where  U  is  the  wind  speed  in  m/s,  9  is  the  local  incidence  angle  in  degrees,  and  cp  is  the  relative 
wind  direction  (defined  as  zero  blowing  towards  the  radar).  The  correction  is  derived  from  clutter 
data  that  is  limited  to  observations  at  two  wind  speeds  (7  and  13  m/s)  over  a  rather  narrow 
incidence  angle  range  (35°  to  50°),  and  is  assumed  to  be  independent  of  wind  aspect  angle. 
However,  the  correction  will  be  used  to  convert  VV  to  cross-polarization  clutter  at  all  wind 
speeds  and  incidence  angles  considered.  Representative  cross-polarization  clutter  transects  based 
upon  this  model  modification  are  also  included  in  Figure  2. 

The  K-distribution  is  used  to  describe  the  distribution  of  digital  numbers  in  the  image  due  to 
image  fading  (speckle)  and  local  modulation  of  the  mean  clutter  level  by  meteorological  or 
oceanographic  processes.  The  K-distribution  includes  two  parameters:  the  ENL  describes  the 
speckle  fading;  the  order  parameter  v  describes  the  local  modulation  of  the  mean.  For  large  order 
parameters,  the  K-distribution  converges  on  the  chi-squared  distribution  (i.e.,  the  clutter  becomes 
Gaussian).  For  small  order  parameters,  the  clutter  becomes  rather  inhomogeneous,  and  the  tail  of 
the  K-distribution  becomes  long. 

A  key  problem  is  what  order  parameter  to  use  for  purposes  of  ship  detectability  modelling.  Some 
empirical  estimates  of  the  order  parameter  from  an  ERS  and  RADARSAT-1  database  of  images 
taken  over  operational  buoys  are  shown  in  Figure  3  [31].  Although  it  is  desirable  to  model  the 
order  parameter  as  a  function  of  the  wind  speed  (e.g.,  the  model  shown  as  a  blue  line  in  the 
figure),  in  general,  this  sort  of  dependence  is  not  observed.  Although  there  appear  to  be  fewer 
occurrences  of  large  order  parameters  as  the  wind  speed  increases,  it  is  seen  that  the  order 
parameter  can  be  small  at  any  wind  speed. 

3.2.3  Ship  RCS  and  its  variability 

There  is  rather  little  discussion  of  ship  RCS  in  the  open  literature.  In  response  to  requirements  to 
characterize  RADARSAT-1’ s  ship  detection  performance,  a  semi-empirical  model  was  developed 
based  upon  a  very  limited  number  of  observations  of  known  ships  [30].  The  initial  model  was 
parameterized  in  terms  of  the  ship  length  and  the  local  incidence  angle.  The  model  applied  to  HH 
polarization  only  and  was  considered  to  provide  a  useful  rule-of-thumb. 

With  the  advent  of  automatic  identification  system  (AIS),  a  mandated  VHF  transponder  system 
for  collision  avoidance  in  commercial  shipping,  it  has  become  feasible  to  generate  very  large 
databases  of  validated  ship  signatures.  At  DRDC  Ottawa,  shore -based  commercial  AIS  data  were 
obtained  in  conjunction  with  several  RADARSAT-1  and  Envisat  ASAR  trials  [32]  [34]  [36],  for 
compilation  of  a  database  of  more  than  4000  validated  ship  signatures  that  may  be  used  to 
improve  models  of  ship  RCS  and  its  variability. 
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R-1 ,  ERS-1/2:  N=131 


Figure  3:  K-distribution  order  parameter  as  a  function  of  wind  speed  for  open  ocean  conditions. 

For  Sentinel- 1  modelling,  the  focus  is  on  Envisat  ASAR  ship  signatures  since  they  provide 
polarization  information.  The  ship  RCS  data  were  evaluated  to  establish  parametric  dependencies 
[36].  The  only  strong  correlation  that  was  noted  was  on  the  ship  length;  no  systematic 
dependence  on  incidence  angle  was  observed. 

Based  upon  the  AIS  database,  the  initial  model  for  RADARSAT-1  data  established  a  very  simple 
relationship  of  the  form: 

cr  =  1.25  L2  (2) 

where  o  is  ship  RCS  in  square  metres  and  L  is  the  ship  length  in  metres.  This  model  was  very 
consistent  with  the  earlier  rule-of-thumb  model,  with  the  variability  of  the  ship  RCS  normalized 
by  equation  (2)  then  describing  the  fading  of  the  ship  RCS  due  to  un-modelled  parameters  such  as 
ship  aspect  angle,  local  incidence  angle,  sea  state,  ship  type,  etc. 

However,  it  was  noticed  that  the  resulting  PDF  did  not  have  unity  area  or  mean.  For  purposes  of 
the  Sentinel- 1  analysis,  the  model  fits  to  the  Envisat  ASAR  data  were  re-visited  such  that  the  PDF 
was  forced  to  have  unity  area  and  mean.  This  requirement  was  constrained  by  retaining  the 
quadratic  ship  length  dependence  that  was  found  previously.  The  generalized  model  is  then 
described  as: 
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(3) 


er pp  A. pp  L 


where  APP  is  a  fitted  coefficient  and  PP  represents  the  polarization  (H  or  V)  on  transmit  and 
receive.  Values  for  APP  are  tabulated  in  Table  5  for  the  various  polarization  choices  and  with 
stratification  of  the  data  by  ENL,  which  is  range  dependent  for  Envisat  ASAR.  In  the  Table,  N  is 
the  number  of  samples  available  for  each  case,  Co-Pol  includes  both  HH  and  VV  data,  and  Cross- 
Pol  includes  both  HV  and  VH  data. 

To  model  the  impact  of  the  choice  of  PD,  the  variability  of  the  normalized  ship  RCS  was  used  to 
establish  the  10th  and  20th  percentile  (lOp  and  20p)  values.  Numerous  attempts  were  made  to 
model  the  observed  variability  with  simple  probability  density  functions  (PDFs)  using  goodness- 
of-fit  tests  [36].  Although  the  log-normal  and  Gamma  distributions  could  fit  the  observed  data,  it 
was  found  that  the  fit  in  the  vicinity  of  the  tails  of  the  distributions  was  variable  and  generally 
poor.  Rather  than  use  a  fitted  PDF,  the  observed  empirical  cumulative  density  function  (ECDF), 
as  plotted  in  Figure  4,  was  used  instead.  The  lOp  and  20p  levels  can  be  read  off  of  these  curves 
directly,  as  tabulated  in  Table  5.  From  the  plots  and  the  table,  it  can  be  noted  that  HH  and  VV 
have  similar  distributions,  that  VH  and  HV  have  similar  distributions,  and  that  the  normalized 
RCS  of  the  co-polarization  channels  is  more  dispersed  than  that  of  the  cross-polarization 
channels. 

As  a  check,  the  dependence  on  ENF  was  investigated.  Such  dependence  seems  possible  since  the 
variability  could  be  a  sampling  issue  related  to  ENF.  The  data  appear  in  Table  5  and  have  been 
plotted  in  Figure  5.  It  could  be  argued  that  there  is  a  slight  dependence  on  ENF,  but  the 
difference  between  co-polarization  and  cross-polarization  percentiles  is  much  larger.  As  such,  it 
was  decided  to  model  the  normalized  ship  RCS  PDF  as  independent  of  ENF  and  for  simplicity,  to 
focus  only  on  the  generic  co-polarization  and  cross-polarization  data,  as  appear  in  the  highlighted 
lines  of  Table  5.  For  simplicity,  only  lOp  and  20p  levels,  corresponding  to  PD  =  0.9  and  PD  =  0.8, 
respectively,  will  be  considered. 

Table  5:  RCS  fits  and  ECDF  percentile  values  for  Envisat  ASAR  data.  (Continued  on  next  page). 


^  ^  ^  ^Q^PqI 

ENL 

N 

A 

lOp 

20p 

ALL 

1330 

2.340 

0.154 

0.231 

1.75 

407 

2.851 

0.099 

0.187 

2.25 

222 

2.009 

0.161 

0.242 

2.66 

291 

1.729 

0.182 

0.289 

3.30 

50 

2.391 

0.141 

0.204 

3.75 

360 

2.451 

0.189 

0.279 
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Table  5:  Concluded. 


***  Cross-Pol  *** 

ENL 

N 

A 

lOp 

20p 

ALL 

783 

0.155 

0.284 

0.432 

1.75 

299 

0.176 

0.229 

0.377 

2.25 

158 

0.145 

0.337 

0.460 

2.66 

124 

0.137 

0.296 

0.455 

3.30 

50 

0.118 

0.334 

0.455 

3.75 

152 

0.154 

0.352 

0.462 

***  HH-Pol  *** 

ENL 

N 

A 

lOp 

20p 

ALL 

672 

2.000 

0.152 

0.239 

1.75 

192 

2.248 

0.108 

0.196 

2.25 

107 

1.803 

0.151 

0.213 

2.66 

159 

1.655 

0.156 

0.270 

3.30 

22 

1.863 

0.121 

0.200 

3.75 

192 

2.164 

0.213 

0.307 

***  VV-Pol  *** 

ENL 

N 

A 

lOp 

20p 

ALL 

658 

2.687 

0.155 

0.235 

1.75 

215 

3.390 

0.108 

0.188 

2.25 

115 

2.200 

0.206 

0.288 

2.66 

132 

1.819 

0.216 

0.297 

3.30 

28 

2.806 

0.172 

0.211 

3.75 

168 

2.780 

0.167 

0.249 

***  HV-Pol  *** 

ENL 

N 

A 

lOp 

20p 

ALL 

387 

0.169 

0.258 

0.414 

1.75 

129 

0.195 

0.180 

0.343 

2.25 

74 

0.153 

0.290 

0.422 

2.66 

76 

0.146 

0.324 

0.473 

3.30 

22 

0.098 

0.314 

0.402 

3.75 

86 

0.181 

0.339 

0.464 

***  VH-Pol  *** 

ENL 

N 

A 

lOp 

20p 

ALL 

396 

0.142 

0.310 

0.452 

1.75 

170 

0.161 

0.272 

0.387 

2.25 

84 

0.137 

0.370 

0.484 

2.66 

48 

0.123 

0.297 

0.431 

3.30 

28 

0.134 

0.378 

0.438 

3.75 

66 

0.119 

0.280 

0.506 
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Figure  4:  ECDFs  as  a  function  of  polarization  for  Envisat  ASAR  data:  a)  all  data;  b)  zoom-in  on 

the  lower  tail. 
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ECDF  10  Percentile  vs  ENL 
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Figure  5:  10th  and  20th  percentile  values  of  the  ECDFs  as  a  function  of  ENL  for  Envisat  ASAR 

data. 


3.2.4  Model  implementation 

The  model  elements  are  implemented  as  follows: 

1 .  NESZ  and  ENL  are  used  to  model  the  noise  floor  and  its  PDF  ( Pn ),  which  leads  to  a  noise- 
only  critical  intensity  Icn  that  is  a  function  of  PFA\ 


cn 

(\-PFA)=  J  P„(x)dx 


(4) 


2.  The  wind  speed,  wind  direction,  incidence  angle,  polarization,  ENL,  and  v  are  used  to  model 
the  ocean  clutter  level  and  its  PDF  ( Pc ),  which  leads  to  a  clutter-only  critical  intensity  Icc  that 
is  a  function  of  PFA: 
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(5) 


(\-PFA)=\Pc(x)dx 

0 


3.  The  clutter  and  noise  are  assumed  to  be  statistically  independent  such  that  the  overall  critical 
intensity  is: 

(6) 

4.  The  critical  intensity  is  scaled  upward  by  dividing  by  the  appropriate  lOp  or  20p  threshold 
level  (rjpD),  which  depend  upon  the  radar  polarization,  and  is  then  converted  to  a  ship  RCS  (cr) 
based  upon  the  radar  resolution  in  azimuth  (pA )  and  ground  range  (pgrY 

(j  =  I_cPaPgR_ 

Vpd 


5.  The  ship  RCS  is  adjusted  upward  by  a  margin  of  M  =  3  dB  to  accommodate  the  possible 
impact  of  unmodelled  factors  that  could  degrade  the  peak  target  response,  hence  ship 
detectability  (see  Section  3.2.5); 

6.  The  ship  RCS  is  then  converted  to  an  equivalent  ship  length  ( L )  as  a  function  of  the 
polarization: 


(8) 


The  critical  intensity  calculation  is  computationally  intensive,  thus  a  3 -dimensional  critical 
intensity  matrix  was  derived  and  stored  in  a  look-up  table  as  a  Matlab®  MAT-file.  The  3- 
dimensions  are  represented  by  the  order-parameter  (v),  ENL,  and  the  significance  level  {SigLev  = 
1  -  Pfa),  all  at  high  enough  representative  resolutions  to  ensure  accurate  estimation  of  the  critical 
intensity  level  via  interpolation.  The  dynamic  ranges  of  these  3  parameters  were  configured  with 
the  intention  of  encompassing  all  meaningful  requirements  that  might  arise,  as  represented  by  the 
following  Matlab®  vectors: 

v  =  [2 : 0.5  : 20  22 : 2 : 50  60 : 10 : 100  125  1000]; 

ENL  =  [1:0.1: 8  8.5: 0.5: 15];  (9) 

SigLev  =  1.0- [1(1 0)"1  3(1  Of2  1(1 0)'2  3(1 0)'3  1(10)  3  ...  3(1 0)'10  1(1 0)'10]. 
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Other  specific  values  may  be  established  through  interpolation  within  the  matrix.  A 
representative  example  is  shown  in  Figure  6.  It  should  be  noted  that  this  look-up  table  applies 
equally  well  to  the  noise-only  (i.e.,  large  v)  and  clutter-only  cases. 

The  model  has  been  implemented  in  Matlab®  as  ship  detectability  analysis  tool,  including  a 
graphical  user  interface,  as  described  in  Annex  A.  All  model  plots  and  tables  presented  in  this 
report  were  generated  using  the  tool. 


Ic  at  SigLevel=1 .0-2.5e-9 


Figure  6:  Visualization  of  I c  for  the  case  ofPFA  =  2.5(1  Of9. 


3.2.5  Model  caveats 

3.2.5. 1  Relative  versus  absolute  model  performance 

Although  the  ship  detection  model  provides  an  estimate  of  the  minimum  detectable  ship  length,  it 
should  be  noted  that  this  is  a  combined  statistical  and  semi-empirical  model,  and  not  all  factors 
have  been  or  can  be  fully  taken  into  consideration.  Since  the  ship  detection  performance  is 
largely  based  upon  the  clutter  level  and  the  SAR  image  quality,  the  model  is  expected  to  be  quite 
adept  at  predicting  the  relative  performance  between  beam  modes  or  radars,  all  else  being  held 
equal.  The  minimum  detectable  ship  lengths  should  be  representative,  but  a  generalized 
validation  is  extremely  difficult,  especially  for  smaller  ships.  Representative  model  predictions  in 
the  context  of  RADARS  AT- 1  and  Envisat  ASAR  data  are  included  in  Annex  B. 
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3.2. 5. 2  Inter-look  coherence 


It  has  been  proposed  that  inter-look  coherence  techniques  could  help  with  the  detection  of  smaller 
ships  [17].  The  premise  is  that  each  “look”,  which  arises  from  bandpass  filtering  of  the  Doppler 
spectrum,  corresponds  to  a  different  observation  time.  In  the  “split-look”  ship  detection 
technique,  two  sub-apertures  are  used.  Ship  targets  of  interest  should  be  correlated  between 
looks,  while  the  background  clutter  remains  uncorrelated.  This  concept  has  been  difficult  to 
verify  and  quantify  [8],  probably  due  to  target  decorrelation  caused  by  target  motion  and  changes 
in  the  target  aspect  angle  between  the  two  looks.  Nevertheless,  a  0.9  dB  improvement  in  contrast 
has  been  noted  for  low  contrast  targets,  i.e.,  signal-to-noise  ratio  (SNR)  <  15  dB  [13]. 

Inter-look  coherence  techniques  or  other  phase-based  analysis  methods  are  not  included  in  the 
ship  detection  model.  To  zeroth  order,  their  impact  could  be  assessed  through  downward 
adjustment  of  the  model  margin.  However,  revision  of  the  clutter  model  parameters  might  also  be 
required. 

3.2. 5.3  Image  pre-processing  to  enhance  ship  targets 

Several  pre-processing  strategies  have  been  proposed  to  enhance  ship  targets  prior  to  detection, 
e.g.,  [27].  The  impacts  of  these  strategies  are  not  included  in  the  ship  detection  model.  To  zeroth 
order,  their  impact  could  be  assessed  through  downward  adjustment  of  the  model  margin. 
However,  revision  of  the  clutter  model  parameters  might  also  be  required. 

3. 2. 5.4  Impulse  response  location  and  target  motion 

There  are  many  other  factors  that  will  impact  the  ship  signature  peak  value  that  have  not  been 
modelled,  including  where  the  peak  value  lies  relative  to  the  sampling  grid  and  the  effects  of 
target  motion,  e.g.,  [21].  To  get  a  handle  on  some  of  these  factors  it  is  assumed  that  the  ship 
signature  is  a  simple  point  target  impulse  response.  The  azimuthal  impulse  response  of  a  moving 
point  target  may  be  written  [29]: 


Sa(  0  = 
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where  t  is  the  azimuth  time,  T  =  Roj3 /V  is  the  coherent  integration  time,  B  =  2V fit/ 1  is  the 
Doppler  bandwidth,  N  is  the  number  of  looks  (i.e.,  ENL),  tD  =  /  F)(v/  E)  scales  the  target 

radial  velocity  component  v  (defined  as  positive  for  target  motion  away  from  the  radar), 
fj  =  Roa  IV2-2uIV  scales  the  target  cross-range  velocity  component  u  and  the  target  radial 
acceleration  component  a  (which  is  also  positive  away  from  the  radar),  /3  =  Al  2pA  is  the 

azimuth  beamwidth,  V  is  the  platform  velocity,  X  is  the  radar  wavelength,  R0  is  the  slant  range  of 
closest  approach,  pA  is  the  1-look  azimuth  resolution,  Gaussian-shaped  functions  have  been  used 
to  model  the  azimuth  antenna  pattern  and  the  look  extraction  filters,  and  linear  FM  and  matched 
filter  processing  is  assumed.  In  equation  (10),  the  first  term  is  the  peak  value  of  the  azimuthal 
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impulse  response,  the  second  term  represents  the  effect  of  Doppler  suppression,  and  the  third  term 
sets  the  shape  and  location  of  the  azimuth  impulse  response.  Strictly  speaking,  this  formula  is 
only  correct  for  the  flat  Earth  case  since  it  contains  the  assumption  that  the  spacecraft  and 
footprint  velocities  are  equal  (the  footprint  velocity  is  usually  10%  smaller  than  the  spacecraft 
velocity),  but  it  is  useful  since  it  allows  certain  parametric  dependencies  to  be  deduced  by 
straightforward  calculation  or  inspection. 

3.2. 5.4.1  Sampling  grid 

An  actual  SAR  system  samples  the  azimuthal  impulse  response  at  the  pulse  repetition  frequency 
(PRF)  fp  >  2 B  ,  where  the  lower  limit  arises  from  the  sampling  theorem.  For  a  static  target  with 

N=  1 ,  the  worst  case  maximum-sampled  value  due  to  the  location  of  the  sampling  grid  relative 
the  location  of  impulse  response  peak  value  occurs  at  a  temporal  shift  of  one-half  of  a  sample  for 
the  lowest  sampling  rate,  i.e.,  at  t  =  1/  4 B  .  This  corresponds  to  a  worst-case  reduction  in  the 
peak  sampled  value  of  -0.9  dB.  This  value  would  change  for  different  azimuth  antenna  and 
matched  filter  patterns  and  if  there  is  a  margin  built  into  the  sampling  design  (i.e.,  the  greater-than 
condition),  but  provides  a  guideline  for  the  impact  of  changing  the  sampling  grid  location. 

3.2. 5.4.2  Doppler  suppression 

Doppler  suppression  arises  when  the  Doppler  spectrum  of  a  moving  target  becomes  shifted 
outside  of  the  Doppler  passband  of  the  processor.  The  processor  is  generally  designed  under  the 
assumption  of  a  static  scene.  By  inspection  of  equation  (1)  and  for  N  =  1,  the  amount  by  which 
the  peak  of  the  impulse  response  is  reduced  is  given  by: 


exp 


(  t2 ) 

_  1 D 

f  v2  4 p2) 

T 2 

V  1  ) 

=  exp 

V1  z2  ) 

(11) 


which  is  plotted  for  nominal  IW  parameters  (i.e.,  V  =  7.5  km/s,  X  =  0.056  m,  pA  =  19.5  m)  in 
Figure  7. 

The  radial  velocity  component  could  arise  from  rectilinear  target  motion  in  the  range  direction. 
To  get  a  handle  on  the  scale  of  the  radial  velocity  component,  consider  a  ship  steaming  at  18 
knots  (9.26  m/s).  The  ship’s  speed  will  be  projected  onto  the  range  direction  through  the  ship 
aspect  angle.  As  a  worst  case,  the  ship  will  be  moving  in  the  range  direction,  either  directly 
towards  or  away  from  the  radar.  In  this  case,  the  rangeward  ship  speed  would  then  be  projected 
onto  the  radial  line  of  sight  through  the  Sine  of  the  local  incidence  angle.  For  Sentinel- 1  IW 
mode,  the  incidence  angle  ranges  from  31°  to  46°,  corresponding  to  radial  velocity  components  of 
4.8  m/s  and  6.8  m/s,  respectively,  with  6.7  m/s  as  the  worst  case  within  the  accessible  IW  swath. 
The  worst  case  is  reduced  if  the  ship  is  not  range  travelling,  and  scales  proportionately  for  other 
ship  operating  speeds.  For  example,  for  the  same  ship  steaming  at  24  knots  (12.35  m/s),  the  worst 
case  radial  velocity  component  would  be  8.9  m/s. 

From  the  plot,  it  is  apparent  that  for  a  quickly  moving,  range-travelling  ship,  Doppler  suppression 
could  reduce  the  peak  response  by  5  dB  or  more. 


20 


DRDC  Ottawa  ECR  2008-218 


Doppler  Suppression  for  IW 


Figure  7:  Doppler  suppression  of  the  azimuth  impulse  response  peak. 

3. 2. 5.4.3  Radial  acceleration  and  cross-range  velocity  suppression 

The  impulse  response  and  in  turn  the  peak  impulse  response  value  will  be  suppressed  by  target 
radial  acceleration  and  target  cross-range  velocity  components.  By  inspection  of  equation  (10) 
and  for  TV  =  1,  the  amount  by  which  the  peak  of  the  impulse  response  is  reduced  due  to  these 
target  motions  is  given  by: 
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The  cross-range  velocity  component  could  arise  from  rectilinear  target  motion  in  the  azimuthal 
direction.  The  effect  is  plotted  for  IW2  mid  swath  parameters  (low-orbit  case)  in  Figure  8. 
Cross-range  velocity  suppression  only  reduces  the  peak  response  by  a  few  hundredths  of  a  decibel 
for  most  practical  ship  speeds,  so  may  be  safely  neglected. 
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Cross-Range  Velocity  Suppression  for  IW2  (mid-swath,  low-orbit) 


Figure  8:  Cross-range  velocity  suppression  of  the  azimuth  impulse  response  peak. 

A  radial  acceleration  component  could  be  imparted  upon  a  ship  that  is  advected  by  an  underlying 
ocean  wave  system.  Smaller  ships  would  be  more  susceptible  to  this  type  of  acceleration  motion. 
As  a  worst  case,  consider  a  Sinusoidal  wave  that  is  travelling  in  the  range  direction.  Then,  it  can 
be  assumed  that  the  maximum  wave  orbital  acceleration  is  imparted  upon  the  ship,  amounting  to: 

a  =  co2A  (13) 


where  CD  —  ->Jgk  (i.e.,  the  deep  water  ocean  wave  dispersion  relationship),  g  is  the  acceleration 
due  to  gravity,  k  =  2n  /  A  is  the  wavenumber,  A  is  the  wave  length,  and  A  is  the  wave  amplitude 

corresponding  to  a  significant  wave  height  ( Hs )  of  2  A2 .  In  practice,  the  radial  component  of  the 
orbital  acceleration  would  be  reduced  if  the  wave  is  other  than  range  travelling.  On  the  other 
hand,  the  acceleration  could  be  increased  if  the  ship  is  moving  into  the  wave  (i.e.,  increasing  the 
encounter  wave  frequency).  Also,  depending  upon  the  transfer  function  between  the  wave  motion 
and  the  ship  motion,  it  is  likely  that  the  orbital  acceleration  imparted  on  the  ship  will  be  smaller 
than  that  suggested  by  the  worst  case  of  equation  (13),  which  is  a  function  of  where  the  ship  lies 
along  the  wave  during  imaging.  An  inherent  assumption  in  this  analysis  is  that  L  «  A ,  but  for 
practical  purposes  it  may  be  assumed  that  L  <  A  /  4  . 
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From  equations  (12)  and  (13),  the  radial  acceleration  peak  suppression  may  be  written: 
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where  (A/  A)  is  a  measure  of  the  wave  slope,  which  in  practical  terms  is  limited  to  1/24  (e.g., 

[19]).  For  IW2  parameters  at  mid  swath  (low  orbit  case),  the  worst  case  radial  acceleration 
suppression  is  plotted  in  Figure  9.  It  is  apparent  that  radial  acceleration  suppression  could  reduce 
the  peak  response  by  several  decibels  for  shorter,  higher  frequency  wave  conditions. 


Radial  Acceleration  Suppression  for  IW2  (mid-swath,  low-orbit) 


Figure  9:  Radial  acceleration  suppression  for  a  wave  of  slope  (A/A). 

To  better  illustrate  the  effect,  several  cases  derived  from  Figure  9  are  presented  in  Table  6.  It  is 
apparent  that  radial  acceleration  suppression  could  cause  a  reduction  of  the  impulse  response  peak 
value  by  several  decibels. 

A  more  complete  analysis  of  the  impact  of  wave-induced  ship  motion  on  SAR  image  focus  would 
require  an  engineering  analysis  of  the  motion  imparted  on  a  ship  by  ocean  waves,  use  of  a 
realistic  wave  spectrum,  and  integration  of  the  target  motion  during  formation  of  the  azimuth 
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aperture.  Accelerating  targets  have  been  discussed  in  the  literature  for  the  case  of  generic 
acceleration  terms  [25]  and  for  heave-induced  ship  motion  [10],  with  the  latter  study  showing 
improvements  in  peak  response  of  3  to  5  dB  if  target  heave  motion  is  taken  into  consideration. 
This  is  consistent  with  the  scale  of  the  radial  acceleration  peak  suppression  estimated  here. 
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[m] 
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(A/A) 
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[m] 
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[dB] 
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100 

0.01 

<25 

-0.74 
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200 

0.05 

<50 

-0.21 
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8 

100 

0.02 

<25 

-2.09 

2 

8 

200 

0.01 

<50 

-0.74 

3 

18 

100 

0.03 

<25 

-3.34 

Table  6:  Worst  case  radial  acceleration  peak  suppression  for  several  sea  states  for  IW2  at  mid 

swath. 


3.2. 5.4.4  Range  smearing 

So  far  only  target  motion  impacts  on  the  azimuthal  impulse  response  have  been  considered. 
However,  a  target  that  is  moving  rapidly  in  the  radial  direction  could  pass  through  several  slant- 
range  resolution  cells  during  the  coherent  integration  time,  causing  the  target  to  be  smeared  in  the 
range  direction,  thus  reducing  the  peak  value  of  the  impulse  response.  This  effect  is  important  if: 
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where  pSR  is  the  slant-range  resolution.  For  IW2  mid  swath  parameters  (low  orbit  case)  and  v  = 
10  m/s,  it  is  apparent  that  the  maximum  value  of  the  left-hand  side  of  equation  (15)  is  around  1.7 
m,  in  comparison  to  pSR  =  2.8  m.  Therefore,  range  smearing  should  not  significantly  degrade  the 
peak  target  value. 

3.2. 5. 5  Ship  signature  size  and  spatial  resolution 

If  the  ship  signature  extends  over  several  resolution  cells,  then  the  peak  contrast  could  be  reduced, 
which  could  in  turn  reduce  vessel  detectability.  In  principle,  the  best  detection  performance 
should  arise  if  the  resolution  cell  size  matches  the  scale  of  the  vessel  that  is  being  detected.  The 
ship  detection  model  assumes  that  the  ship  is  restricted  to  a  single  resolution  cell  in  spatial  extent. 
This  will  rarely  be  the  case,  especially  for  radar  such  as  that  proposed  for  Sentinel- 1,  which  has  a 
rather  high  spatial  resolution.  However,  many  ship  signatures  are  characterized  by  a  single 
dominant  scatterer.  As  such,  the  relative  sizes  of  the  ship  signature  and  the  resolution  cell  are  not 
expected  to  be  a  limitation  of  the  model. 
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Range  Smearing  for  IW2  (mid-swath  low-orbit) 


radial  velocity  v  [m/s] 


Figure  10:  Range  smearing  due  to  target  radial  velocity. 

3. 2. 5. 6  False  alarm  rate  reduction 

Many  practical  ship  detection  systems  employ  false  alarm  rate  reduction  techniques  to  reject 
certain  candidate  ship  signatures  from  further  consideration.  These  techniques  can  include:  target 
clustering,  in  which  isolated  targets  are  joined  into  a  single  target  using  target  scale  restrictions; 
estimation  of  the  ship  length  for  signatures  that  span  multiple  resolution  cells;  estimation  of  the 
target  RCS  in  connection  with  the  ship  signature  length  (e.g.,  through  models  such  as  those  in 
equations  (2)  or  (3));  searching  for  a  ship  wake  signature  in  the  vicinity  of  the  candidate  ship 
signature  as  confirmation  of  a  ship  target;  searching  for  the  presence  of  azimuth  and  range 
ambiguities  for  each  candidate  ship  target  and  rejecting  target  ambiguities  that  presented  as 
candidate  ship  signatures;  and  validating  candidate  ship  signatures  using  external  operational  ship 
information  such  as  AIS  data. 

These  false  alarm  rate  reduction  techniques  could  in  effect  reduce  the  false  alarm,  or  increase  the 
probability  of  missed  detection.  The  impact  of  these  false  alarm  rate  reduction  techniques  has  not 
been  included  in  the  ship  detection  model. 
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3.3  Minimum  detectable  ship  length 

The  minimum  detectable  ship  length  calculations  are  based  upon  the  following  fixed  parameter 
choices: 

•  \J=  12  m/s; 

•  q>  =  0°; 

•  v  =  4; 

•  Pfa  =  2.5(10)'9; 

•  Pd  =  0.9; 

•  M=  3  dB; 

•  Co-pol  and  Cross-pol  operation,  with  HH  considered  for  co-pol. 

The  chosen  wind  speed  is  a  representative  case  that  corresponds  roughly  to  sea  state  5,  and  is  an 
operational  upper  limit  in  which  ship  detection  is  desirable.  A  lower  wind  speed  would  result  in 
better  ship  detection  performance. 

The  chosen  relative  wind  direction  is  for  the  case  of  wind  blowing  towards  the  radar.  Towards 
and  away  from  the  radar  are  the  worst  clutter  cases.  Choice  of  a  different  wind  direction  would 
result  in  better  ship  detection  performance. 

The  chosen  order  parameter  is  a  fairly  pessimistic  case,  but  is  a  value  that  arises  at  all  wind 
speeds.  A  larger  order  parameter  would  result  in  better  ship  detection  performance. 

The  chosen  probability  of  false  alarm  is  rather  low  and  represents  less  than  one  false  alarm  for  a 
500  km  by  500  km  scene.  A  higher  false  alarm  rate  permits  the  detection  of  smaller  ships. 

The  chosen  probability  of  detection  represents  the  performance  of  a  useful  ship  detection  system. 
A  lower  probability  of  detection  precludes  the  detection  of  smaller  ships. 

Some  swath-oriented  results  from  the  minimum  detectable  ship  length  calculations  are  shown  in 
Figure  1 1  through  Figure  14.  The  fixed  parameters  used  in  each  case  are  shown  across  the  title  of 
each  plot.  The  plots  use  the  actual  ground  range  resolution  and  NESZ  specific  to  each  position  in 
the  swath.  The  plots  are  colour  coded  by  sensor  and  include  Sentinel- 1  (blue),  RADARSAT-1 
(red),  and  Envisat  ASAR  (green).  The  beams  are  plotted  alternately  with  solid  and  dashed  lines 
for  each  sensor  and  the  beam  mode  labels  according  to  Table  1  through  Table  4  are  placed  at  the 
centre  of  each  swath  (some  of  the  plots  are  quite  busy).  The  same  Sentinel- 1  data  results  are 
plotted  for  the  HH  and  HV  cases  for  reference  purposes.  The  curves  for  the  HH  case  tend  to  be 
smooth  and  co-linear  due  to  the  clutter-limited  imaging.  The  curves  for  the  HV  case  follow  the 
shape  of  the  noise  floor  due  to  the  noise  limited  imaging.  A  few  cases  have  been  tabulated  at  mid 
swath  (39.1°  incidence  angle,  low-orbit  case)  to  permit  a  more  detailed  performance  comparison 
(see  Table  7). 

By  inspection  of  the  table  and  the  plots,  the  following  is  noted: 
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1.  The  minimum  detectable  ship  length  for  the  groups  of  single  beam  modes  decreases  with 
increasing  incidence  angle,  which  illustrates  the  importance  of  incidence  angle  in  reducing 
the  background  ocean  clutter  level. 

2.  The  minimum  detectable  ship  length  decreases  for  modes  with  increasing  (i.e.,  smaller) 
resolution. 

3.  Comparison  of  the  co-polarization  and  cross-polarization  plots  shows  the  impact  of  changing 
the  polarization.  For  co-polarization,  the  ship  RCS  is  larger,  but  the  clutter  level  is  higher, 
especially  for  smaller  incidence  angles.  The  co-polarization  case  is  generally  clutter  limited. 
For  cross-polarization  the  ship  RCS  is  smaller,  but  the  clutter  is  near  or  below  the  noise  floor; 
the  cross-polarization  case  is  generally  noise  limited.  There  are  significant  benefits  to  using 
cross-polarization,  especially  for  acquisitions  at  smaller  incidence  angles.  Furthermore, 
cross-polarization  provides  more  uniform  ship  detection  performance  across  the  image  swath 
since  the  noise  floor  is  somewhat  independent  of  incidence  angle. 

4.  Comparison  of  the  co-polarization  data  for  PD  =  0.8  and  PD  =  0.9  shows  how  an  increase  in 
this  parameter  decreases  the  detection  performance  for  smaller  ships. 

5.  Considering  Sentinel- 1,  due  to  the  higher  resolution,  there  is  significantly  better  ship 
detection  performance  for  IW  than  RADARSAT’s  SCNB,  a  mode  that  is  being  used  for 
operational  ship  detection,  over  a  comparable  swath  width.  We  also  see  that  there  is  a  penalty 
in  using  EW  for  its  broader  swath  due  to  EW  having  coarser  spatial  resolution  than  IW. 

3.4  Model  sensitivity  analysis 

A  sensitivity  analysis  was  carried  out  to  establish  the  impact  of  the  input  parameters.  This  is 
particularly  useful  since  some  of  the  inputs  are  not  well  known  (e.g.,  the  K-distribution  order 
parameter).  The  baseline  used  for  this  analysis  is  IW2  at  mid  swath  (low  orbit  case)  with  HH 
polarization  and  PD  =  0.9,  along  with  the  other  nominal  parameters  that  were  used  in  the  previous 
section.  Perturbations  to  those  parameters  were  used  to  arrive  at  a  revised  minimum  detectable 
ship  length  and  its  local  rate  of  change  (if  applicable).  Figure  15  and  Table  8  through  Table  16 
summarize  the  analysis  results.  In  the  tables,  the  highlighted  rows  indicate  cases  with  the 
nominal  parameters  that  were  used  in  the  model  runs  of  the  previous  section. 

Figure  15a)  and  Table  8  provide  the  model  sensitivity  to  wind  speed,  showing  how  smaller  ships 
may  be  detected  at  smaller  wind  speeds,  and  that  the  wind  speed  sensitivity  increases  with 
increasing  wind  speed.  The  wind  drives  the  background  clutter  level,  which  increases  with 
increasing  wind  speed. 

Figure  15b)  and  Table  9  provide  the  model  sensitivity  to  wind  direction,  showing  how  smaller 
ships  may  be  detected  for  wind  directions  other  than  towards  (0°)  or  away  from  (180°)  the  radar, 
with  the  greatest  wind  direction  sensitivity  at  angles  that  lie  between  up-wind  and  cross-wind. 
The  wind  direction  drives  the  background  clutter  level,  which  is  maximum  for  the  wind  blowing 
towards  or  away  from  the  radar,  and  is  minimum  for  the  wind  blowing  across  the  radar  look 
direction. 
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HH,  L/=1 2ms-1 ;  (t>=0°;  v=4;  PFA=2.5e-09,  PD=0.9,  Margin=3dB 


Figure  11:  Minimum  detectable  ship  length  for  Sentinel- 1  low  orbit  HH  (blue)  and  RADARSAT-1 

(red). 


HH,  L/=1 2ms-1 ;  (t>=0°;  v=4;  PFA=2.5e-09,  PD=0.9,  Margin=3dB 


Figure  12:  Minimum  detectable  ship  length  for  Sentinel- 1  low  orbit  HH  (blue)  and  Envisat  ASAR 

HH  (green). 
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HV,  L/=12ms-1;4)=0o;v=4;  PFA=2.5e-09,  PD=0.9,  Margin=3dB 


Figure  13:  Minimum  detectable  ship  length  for  Sentinel-1  low  orbit  HV  (blue)  and  RADARSAT-1 

(red). 


HV,  L/=1 2ms-1 ;  (t>=0°;  v=4;  PFA=2.5e-09,  PD=0.9,  Margin=3dB 


Figure  14:  Minimum  detectable  ship  length  for  Sentinel- 1  low  orbit  HV  (blue)  and  Envisat  ASAR 

HV  (green). 
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Table  7:  Minimum  detectable  ship  length  at  mid  swath  for  several  scenarios. 


Sensor 

Mode 

HH 

PD  =  0.8 

M 

HV 

PD  =  0.8 

M 

HH 

PD  =  0.9 
[m] 

HV 

PD  =  0.9 
[m] 

RADARSAT-1 

SI 

149.5 

N/A 

183.1 

N/A 

S4 

54.4 

N/A 

66.7 

N/A 

S7 

30.1 

N/A 

36.8 

N/A 

RADARSAT-1 

FI 

24.1 

N/A 

29.5 

N/A 

F3 

18.5 

N/A 

22.7 

N/A 

F5 

15.5 

N/A 

18.9 

N/A 

RADARSAT-1 

SCNB-near  edge 

197.0 

N/A 

241.3 

N/A 

SCNB-far  edge 

71.7 

N/A 

87.8 

N/A 

Envisat  ASAR 

IS1 

309.1 

132.0 

378.6 

161.6 

IS4 

73.1 

63.8 

89.5 

78.1 

IS7 

39.0 

47.1 

47.8 

57.7 

Sentinel- 1 

SI 

38.0 

18.6 

46.5 

22.8 

S3 

19.3 

13.6 

23.6 

16.7 

S6 

10.0 

11.2 

12.3 

13.8 

Sentinel- 1 

IW1 

34.7 

28.0 

42.5 

34.2 

IW2 

25.6 

24.7 

31.4 

30.3 

IW3 

20.6 

22.7 

25.3 

27.8 

Sentinel- 1 

EW1 

180.9 

95.9 

221.6 

117.4 

EW3 

80.7 

64.7 

98.9 

79.3 

EW5 

55.8 

61.2 

68.3 

74.9 

Figure  15c)  and  Table  10  provide  the  model  sensitivity  to  the  margin  parameter,  showing  how  the 
minimum  detectable  ship  length  scales  with  the  chosen  margin  value.  The  margin  value  could  be 
significant,  depending  upon  how  the  factors  discussed  in  Section  3.2.5  accumulate.  It  could  be 
appropriate  to  increase  the  margin,  especially  if  the  interest  is  in  the  detection  performance  for 
smaller  ships. 

Figure  15d)  and  Table  11  provide  the  model  sensitivity  to  the  K-distribution  order  parameter, 
showing  how  the  minimum  detectable  ship  length  decreases  for  increasing  order  parameter  (i.e., 
for  more  uniform  or  Gaussian  clutter),  and  how  the  sensitivity  increases  for  smaller  order 
parameters. 

Table  12  provides  the  model  sensitivity  to  the  local  incidence  angle,  showing  how  the  minimum 
detectable  ship  length  decreases  for  increasing  incidence  angle.  This  arises  since  the  background 
ocean  clutter  level  decreases  with  increasing  incidence  angle.  This  particular  case  is  clutter 
limited. 
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Table  13  provides  the  model  sensitivity  to  several  different  model  and  system  parameters. 
Considering  the  cases  of  co-polarization,  the  minimum  detectable  ship  length  is  smaller  for  HH 
polarization  than  for  VV  polarization  since  the  HH  ocean  clutter  is  smaller  than  the  VV  ocean 
clutter.  For  cross-polarization,  there  is  only  a  marginal  benefit  in  this  case  since  the  incidence 
angle  is  rather  large.  The  benefit  of  cross-polarization  would  be  more  pronounced  at  smaller 
incidence  angles.  Considering  the  cases  of  the  Sentinel- 1  low  versus  high  orbit,  we  see  a 
marginal  benefit  for  the  low  orbit  case,  arising  mainly  from  the  slightly  larger  incidence  angle. 
Considering  the  cases  of  probability  of  detection,  the  choice  of  0.8  permits  the  detection  of 
smaller  ships,  but  at  the  expense  of  more  missed  detections. 

Figure  15e)  and  Table  14  provide  the  model  sensitivity  to  the  probability  of  false  alarm.  Smaller 
ships  may  be  detected  as  the  false  alarm  rate  is  relaxed,  but  at  the  expense  of  more  false 
detections. 

Table  15  provides  the  model  sensitivity  to  some  of  the  basic  sensor  parameters.  The  model  is 
more  sensitive  to  the  range  resolution  than  the  azimuth  resolution  due  to  their  relative  scales 
(azimuth  resolution  is  much  larger  than  range  resolution).  The  model  is  quite  sensitive  to  ENL, 
but  it  is  only  realistic  to  consider  this  sensitivity  if  ENL  is  traded-off  against  one  of  the  resolution 
values.  In  this  case,  the  model  is  not  sensitive  to  NESZ  since  this  is  a  clutter-limited  scenario. 

Table  16  considers  the  impact  of  a  trade-off  between  ENL  and  the  slant-range  resolution. 
Comparing  the  two  single  look  cases,  the  slant-range  resolution  is  doubled  in  the  second  case, 
representing  the  ship  detection  capability  of  two  individual  looks,  or  perhaps  two  views  of  the 
scene  as  could  arise  from  successive  TOPSAR  [3]  scans,  the  scanning  methodology  that  will  be 
employed  for  Sentinel- 1.  The  idea  here  is  to  trade-off  spatial  resolution  for  two  separate, 
temporally  offset  views  of  the  scene  that  could  be  used  to  estimate  the  ship  velocity.  In  this  case, 
the  minimum  detectable  ship  length  in  each  look  degrades  from  31.4  m  to  44.4  m.  This  could  be 
a  useful  trade-off  if  the  ship  velocity  could  then  be  estimated.  The  high  spatial  resolution  of 
Sentinel- 1  makes  this  trade-off  feasible.  The  subsequent  cases  in  the  table  show  the  sequential 
trade-off  as  the  number  of  looks  increases  and  the  spatial  resolution  decreases  (i.e.,  the 
information  content  is  conserved).  A  degradation  of  minimum  detectable  ship  length  is  apparent 
with  increasing  ENL. 
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SENTINEL-1  ,IW2L,HH,6.nc=39.1,pa=19.5,p  =4.426, NESZ=-27.41,NLKS=1,Lfc12m/s,v=4,PFA=2.5e-09,PD=0.9,M=3dB 


Wind  Aspect  Angle  [deg] 


SENTINEL-1, IW2L,HH,e.nc=39.1,pa=19.5,p  =4.426, NESZ=-27.41  ,NLKS=1  ,U=12m/s,<|>=0o,PFA=2.5e-09,PD=0.9,M=3  dB 


K-Distribution  Order  Parameter 


SENTINEL-1, IW2L,HH,e(.nc=39.1,pa=19.5,p  =4.426, NESZ=-27.41,NLKS=1,U=12m/s,(|)=0°,v=4,PD=0.9,M=3  dB 


Figure  15:  Minimum  detectable  ship  length  as  a  function  of:  a)  wind  speed;  b)  wind  direction;  c) 
model  margin;  d)  order  parameter;  and  e)  probability  of  false  alarm. 
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Table  8:  Ship  detectability  -  sensitivity  to  wind  speed  (U). 


u 

L 

dL/dU 

[m/s] 

[m] 

[si 

3 

11.8 

1.92 

6 

17.8 

2.09 

9 

24.3 

2.26 

12 

31.4 

2.44 

15 

39.0 

2.65 

Table  9:  Ship  detectability  -  sensitivity  to  wind  direction  (cp). 


9 

L 

dL!d(p 

[deg] 

[m] 

[m/deg] 

0 

31.4 

-0.003 

45 

25.8 

-0.197 

90 

17.8 

-0.024 

135 

24.0 

0.163 

180 

29.4 

-0.003 

Table  10:  Ship  detectability  -  sensitivity  to  model  margin  (M). 


M 

L 

dL/dM 

[dB] 

[m] 

[m/dB] 

0 

24.3 

2.82 

3 

31.4 

3.98 

6 

48.6 

5.63 

9 

68.6 

7.94 

Table  11:  Ship  detectability  -  sensitivity  to  order  parameter  (v). 


V 

L 

dLldv 

[m] 

[m] 

2 

38.1 

-4.72 

4 

31.4 

-1.77 

10 

25.9 

-0.43 

100 

20.9 

-0.03 

Table  12:  Ship  detectability  -  sensitivity  to  incidence  angle  (0). 


e 

L 

dL/de 

[deg] 

[m] 

[m/deg] 

@Mid  —  2 

35.8 

-2.50 

8 Mid  ~  1 

33.5 

-2.18 

9 Mid 

31.4 

-1.94 

9 Mid  1 

29.6 

-1.70 

Omm  +  2 

28.0 

-1.38 
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Table  13:  Ship  detectability  -  sensitivity  to  polarization  state  (P),  orbit  altitude,  and  probability 

of  detection  (Pd)- 


p 

Orbit 

Pd 

L 

[m] 

HH 

low 

0.8 

25.6 

0.9 

31.4 

high 

0.8 

28.6 

0.9 

35.0 

HV 

low 

0.8 

24.7 

0.9 

30.3 

high 

0.8 

25.7 

0.9 

31.5 

VH 

low 

0.8 

24.7 

0.9 

30.3 

high 

0.8 

25.7 

0.9 

31.5 

VV 

low 

0.8 

35.7 

0.9 

43.7 

high 

0.8 

39.0 

0.9 

47.7 

Table  14:  Ship  detectability  -  sensitivity  to  probability  of  false  alarm  (Pfa)- 


Pfa 

L 

[m] 

2.5(10)'7 

25.7 

2.5(10)-8 

28.6 

2.5(10)'9 

31.4 

2.5(10)'10 

34.4 

Table  15:  Ship  detectability  -  sensitivity  to  sensor  parameters  (if). 


L 

[m] 

dL/dq 

Paz  [m] 

31.4 

0.80  [-] 

Psr  [m] 

31.4 

3.37  [-] 

ENL 

31.4 

-10.91  [m/look] 

NESZ  [dB] 

31.4 

0.07  [m/dBl 
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Table  16:  Ship  detectability  -  trading  Pgr  far  ENL. 


ENL 

Pgr 

L 

M 

[m] 

[m] 

l 

4.4 

31.4 

l 

OO 

OO 

44.4 

2 

OO 

oo 

34.8 

3 

13.2 

37.4 

4 

17.6 

39.7 

5 

22.0 

41.9 

3.5  Implications  for  Sentinel-1 

The  ship  detection  capability  of  the  Sentinel- 1  SAR  was  modelled  by  using  a  statistical  detection 
model  that  includes  system  noise,  ocean  clutter,  and  ship  RCS.  Each  model  element  includes  a 
level  and  its  variability,  which  is  described  by  an  appropriate  PDF.  For  ship  RCS,  a  database  of 
Envisat  ASAR  AP  mode  ship  signatures  that  were  validated  by  shore -based  AIS  data  was  used  to 
derive  a  relationship  between  ship  RCS  and  ship  length,  and  to  model  the  ship  RCS  variability  as 
a  function  of  polarization.  The  model  includes  a  margin  to  represent  imaging  issues  that  could 
not  be  fully  modelled.  The  model  output  is  the  minimum  detectable  ship  length  for  each  of  the 
radars  of  interest,  specifically,  Sentinel- 1,  RADARSAT-1,  RADARS AT-2,  and  Envisat  ASAR. 

The  relative  performance  between  beam  modes,  all  else  held  equal,  is  expected  to  be  a  solid 
prediction  of  the  model.  The  absolute  performance  is  subject  to  certain  assumptions  concerning 
the  environmental  conditions  and  the  radar  performance.  However,  the  absolute  performance  is 
representative  of  the  expected  performance  of  a  particular  radar  or  mode. 

Based  upon  the  analysis  carried  out,  the  following  conclusions  with  respect  to  Sentinel- 1 
capabilities  are  offered: 

•  The  Sentinel- 1  IW  mode  ship  detection  performance  is  significantly  better  than  that  of 
RADARSAT  SCNB  mode  at  HH  polarization  over  a  comparably  wide  swath.  Note  that 
SCNB  is  a  fairly  standard  mode  for  ship  detection  using  RADARSAT  data.  This  improved 
performance  is  due  to  the  much  higher  spatial  resolution  that  will  be  offered  by  IW. 

•  The  Sentinel- 1  IW  mode  ship  detection  performance  at  cross-polarization  is  comparable  to 
the  ship  detection  performance  at  co-polarization,  although  the  cross-polarization 
performance  is  largely  independent  of  incidence  angle.  This  is  because  IW  starts  at  around 
31°  incidence  angle.  The  benefit  of  cross-polarization  is  more  important  for  smaller 
incidence  angles,  as  seen  for  the  performance  of  the  EW  mode. 

•  There  is  a  penalty  in  terms  of  ship  detection  performance  if  the  EW  mode  is  used  to  gain 
additional  swath  coverage.  However,  the  EW  mode  performance  for  the  HV  channel  is 
comparable  to  that  of  RADARSAT  SCNB,  but  with  a  broader  swath. 
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4  Other  maritime  applications 


4.1  Iceberg  detection 

4.1.1  Problem  description 

The  iceberg  detection  problem  is  very  similar  to  the  ship  detection  problem,  except  that  there  is 
often  an  additional  step  taken  to  discriminate  candidate  ship  signatures  from  candidate  iceberg 
signatures,  e.g.,  [10],  [19].  The  discrimination  operation  can  be  based  upon  the  object’s  RCS, 
shape,  and  texture  parameters,  or  on  the  relative  response  of  the  cross-polarization  channel 
compared  to  the  co-polarization  channel. 

Operational  iceberg  detection  tends  to  use  higher-numbered  Envisat  ASAR  AP  mode  imagery. 
Acquisition  in  HH+HV  mode  facilitates  ship/iceberg  discrimination. 

4.1 .2  Implications  for  Sentinel-1 

Sentinel- 1,  with  its  default  selectable  dual  polarization  and  high  spatial  resolution,  will  support 
iceberg  detection  and  ship/iceberg  discrimination.  The  IW  mode  has  a  much  broader  swath  than 
Envisat  ASAR  APP,  which  will  benefit  this  application. 

4.2  Wind  retrieval 

4.2.1  Problem  description 

The  concept  of  kilometre-scale  wind  field  retrieval  from  SAR  ocean  images  was  discussed  in 
Section  3.2.2  in  the  context  of  the  background  ocean  clutter.  Wind  retrieval  from  SAR  is  an 
application  that  has  recently  matured;  SAR-derived  wind  fields  are  being  demonstrated 
operationally  and  are  being  delivered  to  forecaster’s  desktops  in  near-real  time,  providing  positive 
impact  on  weather  forecasts  [14]. 

Aside  from  knowledge  of  the  wind  direction  and  the  SAR  geometry,  a  suitable  polarization- 
dependent  geophysical  model  function  relating  ocean  backscatter  to  wind  speed  is  required.  To 
date,  for  C-band  instruments,  these  tend  to  focus  on  VV  polarization  due  to  model  function 
heritage  to  the  ERS  scatterometers.  Conversion  of  C-band  VV  model  results  to  other 
polarizations  tends  to  be  done  on  an  ad  hoc  basis.  Nevertheless,  errors  of  only  a  few  m/s  have 
been  observed  for  RADARSAT-1  wind  retrieval  for  moderate  wind  speeds  [31].  Cross¬ 
polarization  is  a  bigger  problem  since  there  is  few  validated  observations  and  the  CNR  tends  to  be 
poor;  ad  hoc  HV  or  VH  models  are  immature. 

In  addition  to  an  adequate  CNR,  a  related  issue  is  the  instrument  sensitivity  to  wind  speed.  This 
may  be  assessed  from  the  slope  of  the  geophysical  model  function.  Considering  co-polarization 
only,  the  geophysical  model  function  sensitivity  has  been  plotted  in  Figure  16  for  the  case  of  the 
wind  blowing  towards  the  radar  (cp  =  0°).  Recall  that  the  co-polarization  ratio  is  independent  of 
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wind  speed  and  direction.  Note  that  sensitivity  decreases  for  smaller  incidence  angles  and  for 
higher  wind  speeds.  For  the  latter,  we  note  the  local  maximum  in  this  particular  model  function 
for  higher  wind  speeds. 

Operational  wind  field  retrieval  tends  to  use  broader  swath  ScanSAR  modes  such  as 
RADARSAT-1  SCW  and  Envisat  ASAR  WSM,  to  achieve  wide  area  coverage.  Better 
performance  is  expected  for  smaller  incidence  angles.  HH  or  VV  polarizations  have  been  used 
successfully,  although  VV  polarization  offers  higher  CNR  and  a  more  mature  geophysical  model 
function.  An  impediment  to  operational  utilization  is  infrequent  repeat  coverage  (compared  to  the 
time  scale  of  the  meteorological  processes  of  interest)  and  locally  fixed  pass  times  (i.e.,  always  at 
nearly  the  same  local  time). 


a) 


C-band  Ocean  Clutter:  $=0° 


incidence  angle  0  [deg] 


b) 


C-band  Ocean  Clutter,  Sensitivity  to  Incidence  Angle:  $=0° 


incidence  angle  6  [deg] 


c) 


C-band  Ocean  Clutter:  (|>=0o 


wind  speed  U  [m/s] 


Figure  16:  Wind  retrieval  sensitivity:  a)  clutter  as  a  function  of  incidence  angle;  b)  clutter 
sensitivity  to  incidence  angle;  c)  clutter  as  a  function  of  wind  speed;  d)  clutter  sensitivity  to  wind 

speed. 
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4.2.2  Implications  for  Sentinel-1 


Sentinel- 1  IW  mode  data  appear  to  be  well-suited  to  wind  retrieval  by  using  the  co-polarization 
channel.  There  is  enhanced  sensitivity  to  wind  speed  for  smaller  incidence  angles,  not  to  mention 
the  possibility  of  a  higher  clutter  level,  thus  better  clutter-to-noise  ratio;  an  IW  mode  that  starts  at 
a  smaller  incidence  angle  would  be  beneficial.  If  a  broader  swath  is  needed,  then  the  EW  would 
be  a  good  choice.  However,  the  co-polarization  channel  of  IW  should  be  well  suited  to  wind 
retrieval  in  the  littoral  zone.  A  constellation  of  two  Sentinel- 1  satellites  would  improve  the  repeat 
coverage,  especially  if  taken  in  the  context  of  other  operational  sources  of  SAR  data. 

4.3  Oil  spill  detection 

4.3.1  Problem  description 

It  is  well  known  that  ocean  backscattering  is  a  Bragg  scattering  process  from  ocean  surface 
roughness  at  the  scale  of  the  radar  wavelength,  and  that  a  layer  of  oil  suppresses  the  Bragg-scale 
roughness.  This  in  turn  reduces  the  local  ocean  backscatter,  so  that  oil  can  be  seen  in  a  SAR 
image  as  an  area  of  darker  return.  In  Figure  17,  examples  are  presented  of  backscatter 
suppression  from  two  RADARSAT-1  SCNA  mode  images  that  were  acquired  through  an 
operational  oil  spill  detection  program  in  Canada  [4].  The  lower  portion  of  each  plot  shows  the 
normalized  radar  backscatter  across  the  entire  swath.  The  upper  portion  shows  more  detail  in  the 
vicinity  of  the  oil  slick  of  interest. 

For  the  30  Sept.  2004  image,  there  is  a  backscatter  suppression  of  roughly  5  dB  relative  to  a 
background  clutter  level  of -18  dB.  The  background  clutter  level  corresponds  to  a  retrieved  wind 
speed  of  between  4  m/s  and  8  m/s,  depending  upon  the  wind  direction,  with  the  lower  limit 
corresponding  to  an  upwind  look  direction  and  the  upper  limit  to  a  crosswind  look  direction.  For 
the  23  Nov.  2004  image,  there  is  a  backscatter  suppression  of  roughly  6  dB  relative  to  a 
background  clutter  level  of -17  dB.  The  background  clutter  level  corresponds  to  a  retrieved  wind 
speed  of  between  5  m/s  and  10  m/s,  depending  upon  the  wind  direction. 

In  both  cases,  the  backscatter  has  been  suppressed  to  roughly  -23  dB,  a  few  decibels  above  the 
instrument  noise  floor.  In  these  cases  for  which  the  wind  speeds  were  rather  light,  the  oil  slicks 
were  easy  to  detect  by  visual  inspection  of  the  image  due  to  their  dark  signatures  and  size. 
Detection  relies  both  on  the  reduction  of  the  clutter  and  the  sharpness  of  the  transition  between  oil 
and  the  clean  water  around  the  edge  of  the  spill.  The  images  suggest  that  a  CNR  of  at  least  6  dB 
would  permit  the  reliable  detection  of  oil  slicks.  On  the  other  hand,  if  the  wind  speed  is  too  high, 
the  slick  would  be  dispersed  and  mixed  downward  into  the  water  column.  A  useful  rule  of  thumb 
is  that  the  wind  speed  should  be  no  higher  that  10  m/s  for  successful  oil  spill  detection. 
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Figure  17:  RADARSAT-1  SCNA  C-band  HH polarization  ocean  backscatter  including  an  oil  slick 
off  the  east  coast  of  Canada:  a)  30  Sept.  2004;  b)  23  Nov.  2004. 

Model  results  for  suppression  of  C-band  radar  cross  section  by  a  layer  of  oil  on  the  ocean  surface 
can  help  to  better  justify  the  6  dB  CNR  requirement.  Charbonneau  [2]  has  recently  published 
calculations  of  ocean  backscatter  suppression  that  are  based  upon  modifications  to  an  ocean 
backscatter  model.  The  calculations  combine  an  ocean  surface  wave  density  spectrum  model, 
which  is  used  as  the  “clean”  ocean  reference,  with  the  physical  properties  of  an  oil  layer,  which 
significantly  reduces  the  surface  wave  density  spectrum.  The  amount  of  backscatter  suppression 
is  a  function  of  the  oil’s  physical  properties  (viscosity,  elasticity,  and  surface  tension)  and  of  the 
oil  layer  thickness.  Using  an  oil-damping  model  in  the  clean  surface  spectrum  simulates  the 
reduction  in  the  wave  density  spectrum.  The  oil’s  dielectric  properties  are  included  as  a  classical 
multi-layer  system,  but  have  little  effect  on  the  backscatter  interaction  due  to  the  low  relative 
permittivity  and  the  thickness  of  the  oil  layer  relative  to  the  C-band  wavelength. 

Charbonneau  has  kindly  provided  us  with  several  cases  that  are  relevant  to  this  problem.  These 
are  presented  in  what  follows  as  plots  of  the  backscatter  suppression  as  a  function  of  incidence 
angle  for  several  wind  speeds.  These  backscatter  suppressions  should  be  interpreted  relative  to 
the  background  ocean  clutter  levels  for  a  clean  ocean  surface,  as  are  available  from  the 
appropriate  CMOD  or  hybrid  backscatter  models. 

Figure  18  shows  ocean  backscatter  suppression  for  a  very  thin  (<  10'7  m)  layer  of  Oleyl  Alcohol, 
which  has  similar  properties  to  those  of  a  natural  surfactant.  The  RCS  enhancement  (i.e.,  positive 
degree  of  RCS  suppression)  at  25°  incidence  angle  is  an  artefact  of  modifications  to  the  integral 
equation  model  that  are  meant  to  compensate  for  a  lack  of  signal  near  nadir.  The  prediction  of 
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RCS  enhancement  represents  a  failure  of  this  model  for  this  particular  case;  however,  the  results 
and  trends  are  thought  to  be  correct  for  incidence  angles  that  are  larger  than  35°.  The  increasing 
RCS  suppression  with  increasing  incidence  angle,  and  decreasing  backscatter  suppression  with 
increasing  wind  speed  are  the  key  predictions  of  the  model.  Backscatter  suppressions  of  up  to  7 
dB  are  noted  for  HH  polarization. 

Figure  19  shows  ocean  backscatter  suppression  for  a  20%  emulsion  of  IFO-180,  a  high  viscosity 
fuel  oil.  The  backscatter  enhancement  artefact  at  smaller  incidence  angles  is  no  longer  present, 
although  there  is  a  trend  to  decreasing  backscatter  suppression  with  smaller  incidence  angles. 
There  is  increasing  backscatter  suppression  for  increasing  incidence  angles,  with  backscatter 
suppressions  of  more  than  20  dB  for  lower  wind  speeds,  and  around  10  dB  for  higher  wind 
speeds.  The  backscatter  suppressions  are  similar  for  both  polarizations. 

Figure  20  shows  ocean  backscatter  suppression  for  a  thin  (3xl0-4  m)  layer  of  Iranian  Heavy  crude 
oil,  which  is  medium  viscosity  oil.  We  note  that  the  backscatter  enhancement  artefact  has 
returned  for  VV  polarization.  Otherwise,  we  see  increasing  backscatter  suppression  for 
increasing  incidence  angles,  with  significantly  more  backscatter  suppression  for  HH  polarization 
than  for  VV  polarization. 

These  model  results  affirm  that  the  suggested  minimum  CNR  of  6  dB  is  a  reasonable  requirement 
for  oil  spill  detection.  To  establish  the  impact  of  this  guideline,  consider  the  useful  swath  as  a 
function  of  NESZ.  This  is  plotted  for  Sentinel- 1  IW  mode  parameters  with  a  nominal  NESZ  of 
-24  dB  in  Figure  21.  The  vertical  line  represents  the  upper  limit  wind  speed  of  10  m/s.  It  is 
apparent  that  there  is  more  useful  swath  for  this  application  for  VV  polarization  and  upwind  look 
directions. 

Operational  oil  spill  detection  tends  to  use  RADARSAT-1  SCNA  data  since  much  of  the  swath  is 
suitable  for  oil  detection,  and  higher  resolution  is  available  than  with  RADARSAT-1  SCW  or 
Envisat  ASAR  WSM.  Use  of  VV  polarization  broadens  the  useful  swath;  it  is  expected  that 
operational  services  will  soon  be  based  upon  RADARSAT-2  SCNA  with  VV  polarization  (or 
more  likely  VV+VH  dual  polarization  to  facilitate  ship  detection  in  the  cross-polarization 
channel). 

4.3.2  Implications  for  Sentinel-1 

Although  they  have  comparable  swath  widths,  Sentinel- l’s  IW  mode  provides  a  fairly  limited 
useful  swath  for  oil  spill  detection  compared  to  RADARSAT-1  SCNA  since  it  starts  at  an 
incidence  angle  of  roughly  30°  compared  to  SCNA’s  20°.  Best  performance  will  occur  for  VV 
polarization.  The  useful  swath  for  this  application  would  be  larger  if  IW  started  at  a  smaller 
incidence  angle.  Access  to  two  Sentinel- 1  satellites  would  improve  repeat  coverage.  Availability 
of  dual  polarization  permits  both  ship  detection  and  oil  spill  detection. 
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a) 


Oleyl  Alchohol:  <  1 .0e-007  m,  C-W 


b) 


Oleyl  Alchohol:  <  1 .0e-007  m,  C-HH 


Incidence  Angle  [deg] 


35  40  45  50  55 

Incidence  Angle  [deg] 


Figure  18:  C-band  ocean  backscatter  suppression  by  a  thin  layer  of  Oleyl  Alcohol:  a)  VV;  b)  HH. 


Figure  19:  C-band  ocean  backscatter  suppression  by  IFO180,  20%  emulsion:  a)  VV;  b)  HH. 


Figure  20:  C-band  ocean  backscatter  suppression  by  a  thin  layer  of  Iranian  Heavy  crude  oil:  a) 

VV;  b)  HH. 
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Sentinel-1  Performance  (CNR=6dB,  oNE°=-24db) 


Figure  21:  Swath  having  a  CNR  larger  than  6  dB  for  a  nominal  noise  floor  of -24  dB. 


4.4  Sea  ice  surveillance 

4.4.1  Problem  description 

Use  of  SAR  imagery  for  operational  sea  ice  surveillance  is  a  proven  operational  capability,  e.g., 
[6].  The  radar  backscatter  from  sea  ice  depends  upon  many  complex  and  inter-related  factors 
including  location,  season,  and  ice  type.  Examples  of  backscatter  profiles  from  RADARSAT-1 
SCW  imagery  are  shown  in  Figure  22,  which  were  acquired  near  Anticosti  Island  and  just  north 
of  Newfoundland  during  routine  sea  ice  surveillance  activities.  These  examples  include  the 
smooth  new  ice  type,  which  has  very  low  backscatter.  It  is  apparent  that,  for  this  case,  the 
backscatter  drops  down  towards  and  below  -26  dB,  which  is  the  nominal  noise  floor  for 
RADARSAT-1. 

The  identification  of  smooth  new  ice  is  important  for  smaller  vessels  that  have  limited  capability 
to  operate  within  ice.  Due  to  its  lack  of  mechanical  strength,  smooth  new  ice  deforms  (and 
roughens)  very  quickly  within  the  ice  pack.  As  a  result,  within  dynamic  ice  regimes,  smooth  new 
ice  does  not  exist  over  large  areas  for  extended  periods  of  time.  In  less  dynamic  regions,  such  as 
in  sheltered  bays,  lakes,  and  in  the  lee  of  coastal  areas,  smooth  new  ice  has  a  greater  presence. 
These  sheltered  areas  may  be  high  traffic  areas  for  small  vessels,  thus  it  is  an  important 
requirement  to  be  able  to  reliably  identify  and  map  the  extent  of  smooth  new  ice. 
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a) 


datgl.001,  SCW,  1-FEB-2005  10:19:48.916,  W1W2W3S7,  Descending,  Lstart :  3053,  Lstop  :  3058 


b) 


datgl.001,  SCW,  15-FEB-2005  21:33:28.941,  W1W2W3S7,  Ascending,  Lstart:  3230,  Lstop  :  3246 


Figure  22:  RADARSAT-1  SCW  C-band  HH polarization  sea  ice  backscatter  for  smooth  new  ice: 
a)  near  Anticosti  Island,  1  Feb.  2005;  b)  near  Newfoundland,  15  Feb.  2005. 

To  further  consider  C-band  backscatter  from  sea  ice,  consider  the  multi-polarization  model  of 
Nghiem  and  Bertoia  [15],  which  has  been  validated  for  incidence  angles  up  to  50°  using  JPL 
AirSAR  data.  The  model  was  recently  extended  beyond  50°  and  to  include  the  cross-polarization 
case  [16].  Model  results  are  shown  in  Figure  23;  all  linear  polarizations  are  represented  for  first- 
year  ice  (FY),  multi-year  ice  (MY),  and  smooth  new  ice  (Leads).  It  is  apparent  that  the 
normalized  RCS  difference  between  HH  and  VV  polarization  increases  for  increasing  incidence 
angle,  especially  for  the  case  of  smooth  new  ice.  There  is  good  SNR  for  multi-year  and  first  year 
ice.  There  are  two  potential  problem  cases:  first,  smooth  new  ice  since  the  co-polarization 
backscatter  signature  drops  down  to  and  below  -25  dB,  as  was  noted  in  the  RADARSAT-1  data; 
and  second,  cross-polarization  backscatter,  for  which  aside  from  the  multi-year  ice  case,  the 
backscatter  is  below  -30  dB. 

Ice  edge  detection  is  another  sea  ice  surveillance  requirement.  This  relies  primarily  on  there 
being  a  difference  between  the  sea  ice  backscatter  and  the  ocean  backscatter,  which  is  generally 
most  pronounced  at  larger  incidence  angles.  Of  course,  ice  versus  ocean  texture  also  enters  into 
this,  but  a  reversal  in  ice  and  ocean  backscatter  can  occur,  especially  for  smaller  incidence  angles 
(i.e.,  the  ocean  backscatter  becomes  larger  than  the  ice  backscatter  for  smaller  incidence  angles, 
whereas  the  ocean  backscatter  becomes  smaller  than  the  ice  backscatter  for  larger  incidence 
angles).  Signature  reversal  could  complicate  ice  edge  detection.  Using  the  Nghiem  and  Bertoia 
model,  the  lowest  wind  speed  that  would  cause  a  signature  reversal  (i.e.,  for  the  upwind  look 
direction)  as  a  function  of  the  incidence  angle  at  which  the  signature  reversal  occurs  for  multi- 
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year  ice  and  first  year  ice  is  plotted  in  Figure  24,  which  supports  the  notion  of  signature  reversal 
being  more  likely  for  smaller  incidence  angles. 

Operational  sea  ice  surveillance  tends  to  use  RADARSAT-1  SCW  due  to  its  large  swath.  Dual 
polarization  is  expected  to  be  beneficial  in  some  cases,  and  could  assist  with  ice  edge  detection 
for  smaller  incidence  angles. 


a) 


Ice  Backscatter  (Nghiem  2007):  C-HV 


Figure  23:  C-band  sea  ice  backscatter:  a)  VV;  b)  HH;  c)  HV. 
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Ice/Ocean  signature  cross  over  wind  speed:  Upwind  look  direction 


Figure  24:  Lowest  wind  speed  required  to  induce  ice/ocean  backscatter  signature  reversal. 

4.4.2  Implications  for  Sentinel-1 

In  general,  IW  operating  in  dual  polarization  will  be  well-suited  to  sea  ice  surveillance.  Since 
RADARSAT-1  SCW  mode  is  often  used  for  operational  sea  ice  surveillance,  the  IW  mode  swath 
is  too  narrow  compared  to  standard  operational  practice.  However,  access  to  a  constellation  of 
two  Sentinel- 1  satellites  operating  in  this  mode  would  alleviate  this  issue,  as  could  use  of  EW 
mode. 

The  ability  to  detect  and  classify  the  smooth  new  ice  (or  lead)  class  is  the  most  challenging 
requirement.  Sentinel-l’s  noise  floor  is  similar  to  that  of  RADARSAT-1  or  Envisat  ASAR. 
Also,  the  information  content  of  the  image  (in  terms  of  the  higher  spatial  resolution)  will  be  a 
benefit.  Sentinel- 1  will  perform  similarly  to  RADARSAT-1  or  Envisat  ASAR  in  terms  of  the 
noise  floor,  or  perhaps  better  considering  the  higher  resolution. 

In  Figure  26,  the  vertical  lines  indicate  the  extent  of  the  IW  swath  for  the  low  orbit  case.  For 
multi-year  ice  and  C-band,  HH  polarization,  ice/ocean  backscatter  signature  reversal  can  occur 
for  wind  speeds  in  excess  of  10  m/s.  For  first-year  ice  and  C-band,  HH  polarization,  ice/ocean 
backscatter  signature  reversal  can  occur  for  wind  speeds  in  excess  of  3  m/s.  Signature  reversal 
can  occur  for  first  year  ice,  but  is  less  likely  for  multi-year  ice.  Sentinel- 1  is  less  susceptible  to 
signature  reversal  than  RADARSAT-1  SCW  due  to  the  relatively  large  incidence  angles  that  will 
be  used. 
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4.5 


Ship  wake  detection 


4.5.1  Problem  description 

It  is  well  known  that  ship  wakes  appear  in  SAR  ocean  images  [18],  arising  in  several  forms 
including  classic  Kelvin  and  transverse  wave  wakes  (Figure  25,  e.g.,  [28]),  and  turbulent  ship 
wakes  (Figure  26a).  A  recent  observation  from  Envisat  ASAR  AP  mode  data  is  that  there  can  be 
a  signature  related  to  ship  wakes  in  cross-polarization  SAR  data  (Figure  26b),  perhaps  caused  by 
the  bow  wave  of  the  ship  [9],  [35].  For  the  cross-polarization  case,  the  wake  signature  appears  as 
an  “echo”  of  the  ship  signature  that  is  located  near  the  apex  of  the  wake  that  is  visible  in  the  co¬ 
polarization  channel.  Since  the  echo  is  not  shifted  in  azimuth,  it  cannot  be  part  of  the  ship  itself 
and  must  be  associated  with  the  ship  wake. 

In  general,  ship  wakes  in  co-polarization  imagery  are  most  visible  under  light  to  moderate  wind 
conditions  and  when  the  clutter-to-noise  ratio  is  high,  suggesting  greater  success  for  smaller 
incidence  angles.  While  ship  wakes  could  form  the  basis  of  a  ship  detection  algorithm,  their 
occurrence  is  relatively  rare.  Usually  ship  wakes  are  searched  for  in  the  vicinity  of  candidate  ship 
signatures  as  a  mechanism  to  reduce  false  alarms  and  possibly  to  obtain  an  estimate  of  ship 
velocity  based  upon  the  azimuthal  Doppler  shift  of  the  ship  signature  and  the  wake  orientation. 
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Figure  26:  Envisat  ASAR  image  of  a  ship  and  its  wake  ( ©European  Space  Agency,  2006):  a)  HH 
polarization  turbulent  wake;  b)  HV polarization  wake  signature.  The  arrows  point  to  the  wake 

features  of  interest. 

There  is  much  less  observational  evidence  for  ship  wakes  in  cross-polarization  imagery,  but  they 
do  seem  to  be  more  prevalent  for  smaller  incidence  angles  [9].  In  principle,  for  a  cross¬ 
polarization  wake  signature  it  should  be  possible  to  establish  the  radial  component  of  the  ship 
velocity  from  the  azimuthal  Doppler  shift  of  the  ship  signature.  An  extended  wake  signature  is 
not  apparent  in  the  cross-polarization  case,  so  resolution  of  the  full  velocity  vector  from  a  cross¬ 
polarization  wake  signature  alone  would  be  difficult  unless  the  ship  orientation  can  be  deduced 
from  the  ship  signature. 

Operational  ship  wake  detection  is  rather  rarely  implemented  due  to  opposing  clutter 
requirements  for  ship  wake  detection  and  ship  signature  detection,  with  ship  signature  detection 
taking  precedence  since  ship  wakes  are  rather  rarely  observed.  Dual  polarization  can  help  by 
focussing  ship  detection  on  the  cross-polarization  channel  and  ship  wake  detection  on  the  co¬ 
polarization  channel. 
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4.5.2  Implications  for  Sentinel-1 


Sentinel- 1  in  IW  mode  with  high  spatial  resolution  and  dual  polarization  will  provide  some  ship 
wakes,  especially  if  the  VV  channel  is  available.  Ship  signatures  could  be  sought  in  both 
polarization  channels,  followed  by  wake  detection  in  both  channels  in  the  vicinity  of  the  wake 
signature.  Wake  detection  in  IW  mode  would  be  better  if  smaller  incidence  angles  were 
available. 
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5  Conclusions 


This  document  addresses  the  analysis  of  marine  applications  potential  of  the  GMES  Sentinel- 1 
SAR  mission.  The  main  focus  of  the  work  has  been  on  ship  detectability  modelling,  but  other 
marine  applications  including  iceberg  detection,  wind  retrieval,  oil  spill  detection,  sea  ice 
surveillance,  and  ship  wake  detection  are  also  addressed. 

The  ship  detectability  model  has  extended  previous  modelling  that  was  carried  out  for 
RADARSAT-1  to  other  sensors,  and  takes  into  consideration  alternating  polarization  observations 
from  Envisat  ASAR  trials.  The  detection  model  applies  to  all  linear  polarizations  and  has  been 
implemented  as  a  Matlab®  ship  detectability  performance  tool  that  operates  in  several  modes. 
All  model  results  and  plots  concerning  ship  detection  that  are  presented  in  this  document  were 
generated  with  the  tool. 

Whereas  Sentinel- 1  will  have  several  modes  of  operation,  the  focus  of  this  work  has  been  the 
Interferometric  Wide  (IW)  mode.  This  is  expected  to  be  Sentinel- l’s  default  mode  of  operation 
in  the  littoral  zone.  Performance  of  the  IW  mode  has  been  compared  relative  to  current 
operational  practice  with  RADARSAT-1  and  Envisat  ASAR. 

The  ship  detectability  modelling  has  shown  that  IW  is  very  well-suited  to  ship  detection,  almost 
regardless  of  the  available  polarization.  HH  polarization  is  the  best  of  the  co-polarization 
choices,  with  cross-polarization  improving  performance  at  smaller  incidence  angles.  Availability 
of  two  polarizations  would  improve  the  detection  performance  or  could  be  used  to  reject  false 
alarms.  Cross-polarization  would  have  more  benefits  for  incidence  angles  smaller  than  the 
minimum  of  29°  to  31°  for  IW.  On  the  other  hand,  cross-polarization  clutter  could  become  a 
factor  for  incidence  angles  of  less  than  25°  or  so. 

The  other  maritime  surveillance  applications  were  treated  relatively  cursorily  compared  to  the 
treatment  of  ship  detectability: 

•  For  iceberg  detection  and  ship/iceberg  discrimination,  IW  with  HH+HV  polarization  is 
expected  to  work  well; 

•  For  wind  retrieval,  IW  performance  would  be  better  if  the  incidence  angles  were  smaller; 

•  For  oil  spill  detection,  IW  performance  would  be  better  if  the  incidence  angles  were  smaller; 

•  For  sea  ice  surveillance,  IW  performance  is  expected  to  be  very  good,  but  a  larger  swath 
would  be  desirable; 

•  For  ship  wake  detection,  dual  polarization  with  high  spatial  resolution  is  expected  to  be  a 
benefit,  but  IW  performance  would  be  better  if  the  incidence  angles  were  smaller. 

As  a  general  statement,  the  one-size-fits-all  approach  that  will  be  adopted  with  Sentinel- 1  IW 
operations  will  suit  most  maritime  applications  in  the  littoral  zone  rather  well.  On  the  other 
hand,  use  of  smaller  incidence  angles  would  be  beneficial  for  several  applications. 
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5.1  Future  work  to  consider 


A  few  gaps  and  weaknesses  are  apparent  in  this  study  that  could  benefit  from  additional  study. 

These  are  summarized  as  follows: 

•  Modelling  and  error  bound  estimates  should  be  improved  for  moving  ships; 

•  Cross-polarization  clutter  should  be  measured  and  models  should  be  improved; 

•  Implement  other  C-band  ocean  back  scatter  models  (such  as  CMOD5); 

•  Estimation  of  the  K-distribution  order  parameter  for  ocean  clutter  should  be  addressed  and 
models  for  the  K-distribution  order  parameter  as  a  function  of  sea  state  should  be  developed; 

•  Other  detection  approaches  such  as  use  of  inter-look  coherence  should  be  considered,  both 
experimentally  and  as  an  extension  to  the  current  ship  detectability  tool.  This  would  include 
use  of  successive  TOPSAR  scans  to  provide  two  images  separated  in  time  for  estimation  of 
ship  velocity; 

•  The  RADARS AT-2  component  of  the  ship  detectability  tool  is  incomplete,  especially  for 
the  higher  resolution  RADARSAT-2  modes.  Actual  RADARSAT-2  noise  floors  should 
also  be  implemented  in  the  ship  detectability  tool; 

•  Implement  the  tool  for  other  sensors,  such  as  TerraS AR-X.  This  would  require  a  campaign 
to  acquire  appropriate  data  for  generation  of  ocean  clutter  and  ship  RCS  at  X-band; 

•  RADARSAT-2  UltraFine  mode  data  should  be  evaluated  to  further  study  the  performance 
of  higher  resolution  data,  especially  with  respect  to  the  detection  of  small  ships  (e.g.,  L  <  50 
m)  as  should  be  detectable  by  Sentinel- 1  IW  mode.  UltraFine  mode  data  could  be  degraded 
to  simulate  Sentinel- 1  S  and  IW  mode  data; 

•  Modelling  of  dual-polarization  data  for  ship  and  iceberg  discrimination  should  be  carried 
out  and  implemented  as  an  extension  of  the  ship  detectability  tool. 
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Annex  A  SAR  Ship  Detectability  Performance  Tool 


The  SAR  Ship  Detectability  Performance  Tool  is  a  GUI-driven  system  written  in  the  Matlab® 
scripting  language  and  should  run  on  any  computer  that  has  basic  Matlab®  installed.  The  essence 
of  the  tool  is  to  provide  the  minimum  detectable  ship  length  given  configurable  sensor, 
performance,  and  environmental  parameters.  The  minimum  detectable  ship  length  calculation  is 
described  in  Section  3.2  of  this  document.  The  tool  currently  supports  five  sensors:  Sentinel- 1 
with  Low  and  High  orbit  scenarios,  RADARSAT-1,  RADARS AT-2,  and  Envisat  ASAR. 
Supported  beam  modes  for  each  of  the  sensors  are  shown  in  Table  17.  Other  beam  modes  can  be 
added  and/or  modified  with  relative  ease.  Available  polarizations  are  HH,  VV,  HV,  and  VH.  For 
RADARSAT-1,  only  HH-polarization  is  available,  regardless  of  the  chosen  polarization. 

The  GUI  provides  the  ability  to  modify  parameters  via  a  series  of  pop-up  menus,  radio  buttons, 
and  text  windows.  Text  windows  with  a  white  background  can  be  modified  by  the  user.  Text 
windows  with  a  grey  background  are  not  directly  editable  by  the  user,  but  in  some  cases  are 
slaved  to  another  parameter  that  can  be  modified.  Tool  parameters  are  summarized  in  Table  18. 

Table  1 7:  Supported  beam  modes  for  each  Sensor. 


Sentinel- 1/2 

RADARS  AT- 1/2 

Envisat  ASAR 

SI 

SI 

APP-IS1 

S2 

S2 

APP-IS2 

S3 

S3 

APP-IS3 

S4 

S4 

APP-IS4 

S5 

S5 

APP-IS5 

S6 

S6 

APP-IS6 

IW1 

S7 

APP-IS7 

IW2 

FI 

WSM-SS1 

IW3 

F2 

WSM-SS2 

EW1 

F3 

WSM-SS3 

EW2 

F4 

WSM-SS4 

EW3 

F5 

WSM-SS5 

EW4 

W1 

IMP-IS1 

EW5 

W2 

IMP-IS2 

W3 

IMP-IS3 

EH1 

IMP-IS4 

EH2 

IMP-IS5 

EH3 

IMP-IS6 

EH4 

IMP-IS7 

EH5 

EH6 

ELI 

SCNA 

SCNB 

SCWA 

SCWB 
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Table  18:  Glossary  of  parameters  presented  when  using  the  SAR  ship  detectability  performance 

tool. 


Parameter 

Description 

Status 

Nature 

Units 

Clutter  to  Noise 
Ratio  (CNR) 

The  ratio  of  the  model  ocean  radar  cross- 
section  and  NESZ. 

Output 

dB 

Incidence  Angle 
(0) 

The  angle  between  the  radar  beam  and  the 
local  normal. 

Input 

Sensor 

deg 

K-Distribution 
order  parameter 
(V) 

The  multi-look  intensity  K- distribution 
parameter  that  describes  the  local 
modulation  of  the  mean. 

Input 

Environment 

Margin  ( M) 

A  contingency  factor  that  is  added  to  the 
minimum  detectable  radar  cross  section. 

Input 

Detection 

Performance 

dB 

Minimum 
Detectable  RCS 
(°) 

The  minimum  detectable  radar  cross 
section. 

Output 

dB- 

m2 

Minimum 
Detectable  Ship 
Length  ( L ) 

The  minimum  detectable  ship  length. 

Output 

m 

Minimum 
Detectable  Peak  to 
Clutter  Ratio 
(PCR) 

The  minimum  detectable  peak  to  clutter 
ratio  assuming  the  target  is  restricted  to  a 
single  resolution  cell. 

Output 

dB 

NESZ  (<t°ne) 

The  radar  system  noise  equivalent  sigma 
zero. 

Input  or 
from  Tables 

Sensor 

dB 

NLOOKS  (ENL 
or  N) 

The  equivalent  number  of  independent 
looks. 

Input  or 
from  Tables 

Sensor 

- 

Ocean  Clutter  (<7°) 

The  normalized  radar  cross  section  of  the 
background  ocean  clutter. 

Output 

dB 

PD  %80  or  PD 
%90  (20p  or  lOp) 
(Pd) 

The  two  available  probabilities  of 
detection. 

Input 

Detection 

Performance 

% 

Probability  of 

False  Alarm  ( PFA ) 

The  false  alarm  rate. 

Input 

Detection 

Performance 

- 

RHOa  (pA) 

The  azimuth  resolution. 

Input  or 
from  Tables 

Sensor 

m 

RHOr  ( Pgr ) 

The  ground  range  resolution. 

Input  or 
calculated 

Sensor 

m 

Wind  Speed  (U) 

The  measured  wind  speed,  usually 
parameterized  as  the  equivalent  neutral 
stability  wind  measured  at  a  height  of  10 
m  above  the  ocean  surface. 

Input 

Environmental 

m/s 

Wind  Aspect 

Angle  (<p) 

The  angle  between  the  wind  direction  and 
the  SAR  look  direction  (zero  for  wind 
blowing  towards  the  radar). 

Input 

Environmental 

deg 
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To  install  the  tool,  save  the  contents  of  the  provided  zip  archive  to  a  uniquely  named  directory. 
To  launch  the  tool,  invoke  Matlab®,  change  to  the  installation  directory,  and  type  smspt  at  the 
Matlab®  prompt: 


The  following  startup  menu  is  displayed: 
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The  user  is  presented  with  four  options  from  which  to  choose: 

1.  Sensor-Specific  Ship  Detectability  Calculator; 

2.  Sensor-Specific  Full  Swath  Analysis; 

3.  Generic  Ship  Detectability  Calculator;  and 

4.  Ship  Detectability  Sensitivity  Analysis. 

A.1  Sensor-Specific  Ship  Detectability  Calculator 

Based  on  user-selected  parameters,  the  Sensor-Specific  Ship  Detectability  Calculator  calculates 
the  Minimum  Detectable  Ship  Length  in  meters,  the  Minimum  Detectable  RCS  in  dB-m2,  the 
Clutter  to  Noise  Ratio  in  dB,  the  Ocean  Clutter  level  in  dB,  and  the  Minimum  Peak  to  Clutter 
Ratio  in  dB. 


To  calculate  parameters  for  a  Sentinel- 1  High  Orbit,  HV-polarization,  with  probability  of 
detection  of  80%  for  the  near-edge  of  beam  IW2  at  the  default  supplied  values  of  wind  speed  of 
12  m/s,  wind  aspect  angle  of  0°,  K-distribution  order  parameter  of  4,  probability  of  false  alarm  of 
2.5e-9,  and  margin  of  3  dB,  the  user  would  do  the  following: 


58 


DRDC  Ottawa  ECR  2008-218 


From  the  sensor  popup  menu  select  SENTINEL-1  High  Orbit: 


From  the  polarization  popup  menu  select  HV-Pol: 
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From  the  probability  of  detection  popup  menu  select  PD  80%: 


pt_Se  nso  rS  pecif icCalc  u  lator 


>AR  Ship  Detectability  Performance  Tool 
Sensor-Specific  Calculator 


EL-1  High  Orbit  Q  |HV-Pol 


Select  Beam 


Incidence  Angle  [deg]  M^x  ~  b 

RHOa  [m] 

RHOr  (GrndRng)  [m] 

NESZ  [dB] 

NLOOKS 

warn 

Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 

I  Reset  Defaults 


- Calculated  Results - 

Minimum  Detectable  Ship  Length  [m] 
Minimum  Detectable  RCS  [dB-mA2] 
Clutter  to  Noise  Ratio  [dB] 

Ocean  Clutter  [dB] 

Minimum  Peak  To  Clutter  Ratio  [dB] 


Click  on  the  “Select  Beam”  button  to  initiate  the  beam  selection  popup  menu  and  click-and-drag 
to  IW2\ 


TO- 


SENTINEL-1  High  Orbit 

Select  Beam 


SAR  Ship  Detectability  Performance  Tool 
Sensor-Specific  Calculator 

D  HV-PoT  H  jPD  80"/  “  — - ; - fnmjri 

- M  -  smspt_SelectSentine...  |B>||u||X| 


Incidence  Angle  [deg]  ““;§g| 
RHOa  [m] 

RHOr  (GrndRng)  [m] 

NESZ  [dB] 

NLOOKS 


Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 

I  Reset  Defaults 


Calculate 

iw2w - Calculated  Results - 

g^Hium  Detectable  Ship  Length  [m] 
MjlHium  Detectable  RCS  [dB-mA2] 
r  to  Noise  Ratio  [dB] 

Ocean  Clutter  [dB] 

Minimum  Peak  To  Clutter  Ratio  [dB] 


Click  on  the  “OK”  button  to  return  to  main  page. 
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Next,  click  on  the  Incidence  Angle  text  box  and  change  the  incidence  angle  to  the  near-edge  of 
IW2,  i.e.,  to  34.77°.  The  minimum  and  maximum  incidence  angles  for  the  selected  beam  are 
displayed  just  to  the  left  of  the  Incidence  Angle  text  box  for  informational  purposes.  RHOa, 
RHOr,  NESZ,  and  NLOOKS  are  slaved  to  the  Incidence  Angle  and  their  values  may  change  with 
a  change  to  the  Incidence  Angle. 


AR  Ship  Detectability  Performance  Tool 
Sensor-Specific  Calculator 


By® 


Calculate 


-Calculated  Results- 


Minimum  Detectable  Ship  Length  [m] 
Minimum  Detectable  RCS  [dB-mA2] 
Clutter  to  Noise  Ratio  [dB] 

Ocean  Clutter  [dB] 

Minimum  Peak  To  Clutter  Ratio  [dB] 


Clear 


Finally,  click  on  “Calculate”  to  generate  the  Calculated  Results. 


If  desired,  click  on  the  “Reset  Defaults”  button  to  reset  the  default  values  in  the  lower-left  sector, 
or  click  on  “Clear”  to  clear  the  Calculated  Results.  Note  that  clicking  on  “Calculate”  updates  all 
of  the  calculated  results  (i.e.,  it  is  not  essential  to  click  on  “Clear”  prior  to  “Calculate”). 


DRDC  Ottawa  ECR  2008-218 


61 


A.2  Sensor-Specific  Full  Swath  Analysis 


The  Full  Swath  Analysis  tool  provides  the  user  with  the  ability  to  plot  the  minimum  detectable 
ship  length  ( across  the  full  swath  of  the  beam)  versus  incidence  angle  for  all  the  beams  of  one  or 
more  selected  sensors,  at  the  selected  polarization,  and  at  the  desired  performance  and 
environmental  parameter  settings. 


To  generate  a  plot  of  SENTINEL- 1  Low  Orbit  and  RADARSAT-1  sensors  at  HH-polarization 
with  the  default-supplied  performance  and  environmental  parameters,  click  on  the  RADARSAT-1 
“N”  radio-button  to  change  it  to  “Y” 
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Click  on  the  “Process”  button  to  initiate  processing.  The  processing  status  is  provided  in  the 
Matlab®  workspace: 
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A  plot  of  Minimum  Detectable  Ship  Length  in  meters  (on  a  logarithmic  scale)  versus  Incidence 
Angle  (on  a  linear  scale)  in  degrees  is  produced  (the  title  displays  the  selected  parameters): 


For  these  plots,  the  sensors  are  color  coded  as  follows: 

•  SENTINEL- 1  Low  Orbit  in  Blue  (each  beam  has  an  “L”  qualifier  to  designate  Low  Orbit); 

•  SENTINEL- 1  High  Orbit  in  Blue  (each  beam  has  an  “H”  qualifier  to  designate  High  Orbit); 

•  RADARSAT-1  in  Red  (each  beam  has  an  R1  qualifier  to  designate  RADARSAT-1); 

•  Envisat  ASAR  in  Green;  and 

•  RADARS AT-2  in  Black  (each  beam  has  an  R2  qualifier  to  designate  RADARSAT-2). 

Sequential  beams  are  plotted  as  a  solid  line,  dashed  line,  solid  line,  dashed  line,  etc.  to  help 
differentiate  between  neighbouring  beams.  As  well,  the  beam  mode  name  is  stamped  at  the  mid¬ 
swath  of  each  beam. 

The  “Restore  Defaults”  button  sets  the  sensor  to  SENTINEL- 1,  HH  polarization,  and  the  nominal 
performance  and  environmental  parameters. 
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A.3  Generic  Ship  Detectability  Calculator 


The  Generic  Ship  Detectability  Calculator  is  essentially  the  same  as  the  Sensor  Specific 
Calculator ,  but  without  any  sensor  dependence.  As  such,  all  the  model  parameters  are  editable 
with  no  real  constraints. 


mspt_GenericCaIculator 


SAR  Ship  Detectability  Performance  Tool 
Generic  Calculator 


-  Parameter  Settings- 

jPD  90%  Q 


RHOa  [m] 

RHOr  (GrndRng)  [m] 

Incidence  Angle  [deg] 

NESZ  [dB] 

NLOOKS 
Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 


Reset  Defaults 


Version  1.0 


Calculate 


-Calculated  Results- 


Minimum  Detectable  Ship  Length  [m] 
Minimum  Detectable  RCS  [dB-mA2] 
Clutter  to  Noise  Ratio  [dB] 

Ocean  Clutter  [dB] 

Minimum  Peak  To  Clutter  Ratio  [dB] 


Clear 


At  startup  and  by  clicking  on  the  “Reset  Defaults”  button,  the  parameters  are  set  to  the  nominal 
defaults  for  Sentinel- 1  Low  Orbit  at  mid-swath  of  beam  SI. 

A.4  Ship  Detectability  Sensitivity  Analysis 

The  Sensitivity  Analysis  tool  is  sensor  specific  and  is  designed  to  provide  the  sensitivity  of  the 
minimum  detectable  ship  length  with  respect  to  varying  wind  speed,  wind  aspect  angle,  Re¬ 
distribution  order  parameter,  probability  of  false  alarm,  and  model  Margin. 
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Rssdpt  _ScnsilivilyAiuilyiis 


SAR  Ship  Detectability  Performance  Tool 
Sensitivity  Analysis 


SEHT1HEL-1  Low  Orbit 


|  Select  Beam  | 


Incidence  Angle  [deg] 
RHOa  [m] 

RHOr  (GmdRng)  [m] 
NESZ  [dB] 

NLOOKS 


I4ax  -  2629 
Min- 


Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 


Version  1.0 


Select  Variable  and  Set  LimitS- 

Wtnd  Spend  [m's.]  [f] 

Minimum 
Maximum 
Increment 


Process 


The  user  selects  the  sensor,  polarization,  probability  of  detection  and  the  beam  mode;  the  tool 
automatically  sets  the  incidence  angle,  RHOa,  RHOr,  NESZ  and  NLOOKS  parameters  at  nominal 
mid-swath  settings.  RHOa,  RHOr,  NESZ,  and  NLOOKS  are  still  slaved  to  the  incidence  angle 
but  are  user  editable  if  desired.  The  variable  of  interest  is  selected,  the  dynamic  range  and 
resolution  of  the  selected  variable  are  input  by  the  user,  and  a  plot  of  minimum  detectable  ship 
length  versus  the  selected  variable  is  produced. 

By  way  of  example,  a  wind  aspect  angle  sensitivity  analysis  is  illustrated  for  Sentinel- 1  Low 
Orbit,  HH-pol,  PD  of  90%,  beam  SI  and  nominal  settings  for  wind  speed,  K-distribution  order 
parameter,  probability  of  false  alarm  and  margin.  Prom  the  Select  Variable  popup  menu  select 
the  Wind  Aspect  Angle  parameter: 


»d|it_ScnsititrilyAiui]yxis 


Bi'n&J 


SAR  Ship  Detectability  Performance  Tool 
Sensitivity  Analysis 


|SEHT1MEL-1  Low  Orbit 


Select  Beam 


Incidence  Angle  [deg] 
RHOa  [m] 

RHOr  (GmdRng)  [m] 
NESZ  [dB] 

NLOOKS 


Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 


Select  Variable  and  Set  Limits- 

wind  Spood  [in's]  ff] 


M  Wind  Aspect  Angle  |deg| 

Increment 


Process 
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Next,  set  Minimum ,  Maximum  and  Increment  values  for  the  selected  variable  in  the  editable  text 
boxes  (this  sets  the  variable  resolution  and  dynamic  range  over  which  to  plot  the  minimum 
detectable  ship  length): 


a- 


Esdj>t_Scnsilirilyftiuilysis 


SAR  Ship  Detectability  Performance  Tool 
Sensitivity  Analysis 


SENTINEL-1  Low  Orlilt 


Select  Beam 


Wind  Aspect  Ar 


Maximflm 


Increment 


Incidence  Angle  [deg] 
RHOa  [m] 

RHOr  (GmdRng)  [m] 
NESZ  [dB] 

NLOOKS 


Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 


Select  Variable  and  Set  Limits 


Finally,  click  on  the  “Process”  button  to  produce  a  plot  of  minimum  detectable  ship  length  (on  a 
logarithmic  scale)  versus  wind  aspect  angle  (on  a  linear  scale): 
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All  of  the  parameter  settings  are  provided  in  the  title  and  the  min  and  max  minimum  detectable 
ship  length  in  meters  are  provided  (inset)  for  information  purposes. 

As  another  example,  setting  a  range  of  values  for  the  Probability  of  False  Alarm  can  be  tricky 
since  the  numbers  involved  are  so  small.  However,  choosing  the  parameters  noted  below 
produces  the  following  plot. 


& 


is  d  |i!  _Se  ns  i  I  ivi  I  yi  rud  ysis 


SAR  Ship  Detectability  Performance  Tool 
Sensitivity  Analysis 


SEHT1HEL-1  Low  Orbit 


| Select  Beam] 


Incidence  Angle  [deg] 
RHOa  [m] 

RHOr  (GmdRng)  [m] 
NESZ  [dB] 

NLOOKS 


Max-  26.29 
Min  -  !«■!» 


Wind  Speed  [m/s] 

Wind  Aspect  Angle  [deg] 
K-Distribution  Order  Parameter 
Probability  of  False  Alarm 
Margin  [dB] 


Select  Variable  and  Set  Limits- 

Probability  of  Falso  A...  Q 

Minimum 
Maximum 
Increment 


Process 


^Figure  7 


File  Edit  View  Insert  Tools  Desktop  Window  Help 

D  &  Q  m  k  Q  ^  □  0  Q  Q 


yiDS 


SENTINEL-1  ,S1L  ,HH, 5^=23. 1  ,jO=4.9,p =3.927 ,NESZ=-26.06,NLKS=1  ,U=1 2m/s ,  ^0° ,  ^4  ,PD=0.9 ,M=3  dB 
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Annex  B  Examples  from  Envisat  ASAR  and 
RADARSAT-1 


To  demonstrate  the  relationship  between  the  SAR  Ship  Detectability  Performance  Tool  and  actual 
SAR  data,  several  case  studies  have  been  developed  based  upon  available  Envisat  ASAR  and 
RADARSAT-1  imagery,  with  ship  validation  data  via  coastal  AIS  data.  In  the  Envisat  case,  the 
same  ship  has  been  considered  in  three  different  alternating  polarization  scenes.  In  the 
RADARSAT-1  case,  two  ships  in  a  single  SCNA  scene  have  been  considered.  For  each  case,  the 
overall  image  with  validated  ships  is  presented,  followed  by  target  analysis  using  a  DRDC  Ottawa 
target  analysis  tool  kit.  The  ship  detectability  tool  was  then  run  based  upon  actual  image 
parameters,  and  the  tool  output  was  considered  in  the  context  of  the  available  image  data. 

The  cases  considered  are  summarized  in  Table  19  and  the  tool  outputs  are  provided  in  Table  20. 
In  the  latter  table,  the  wind  direction  was  taken  as  zero,  unless  an  external  wind  direction  source 
was  available.  In  either  case,  the  retrieved  wind  speed  sets  the  CNR  for  the  performance 
calculations.  The  K-distribution  order  parameter  (v)  was  estimated  from  the  actual  number  of 
statistically  independent  looks  for  the  specific  beam  mode  and  the  observed  MSVR  in  the  vicinity 
of  the  ship  of  interest.  The  case  studies  are  presented  and  discussed  in  the  following  sections. 


Table  19:  Demonstration  products  -  case  studies. 


Case 

Date 

Location 

Sensor 

Mode 

Pol 

Ship 

L 

[m] 

e 

[deg] 

PCR 

[dB] 

1 

9  Dec.  2006 

Gibraltar 

ASAR  AP 

ISl 

HH 

Le  Rif, \  Ferry 

129 

15.6 

N/A 

HV 

27.5 

2 

22  Dec.  2006 

Gibraltar 

ASAR  AP 

IS3 

HH 

Le  Rif, \  Ferry 

129 

28.8 

16.6 

HV 

19.9 

3 

16  Dec.  2006 

Gibraltar 

ASAR  AP 

IS7 

HH 

Le  Rif, \  Ferry 

129 

44.4 

30.8 

HV 

23.6 

4 

4  March  2004 

Halifax 

ASAR  AP 

IS3 

HH 

Quest ,  Research 

76 

29.7 

17.7 

VV 

13.1 

HH 

Grenfell,  CCG 

68 

29.8 

16.0 

VV 

12.4 

HH 

Sambro ,  CCG 

16 

29.8 

N/A 

VV 

9.3 

5 

3  Jan.  2007 

Gibraltar 

ASAR  AP 

ISl 

VV 

Nakhchivan ,  Cargo 

140 

14.8 

N/A 

VH 

22.0 

VV 

Stav.  Prince ,  Tanker 

244 

15.0 

15.8 

VH 

32.6 

VV 

Saral ,  Ferry 

119 

15.3 

N/A 

VH 

20.7 

6 

9  July  2006 

Dover 

R1 

SCNA 

HH 

Leo  Glory ,  Tanker 

334 

38.0 

31.3 

HH 

Sandettie ,  Cargo 

142 

36.8 

25.2 
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Table  20:  Demonstration  products  -  ship  detectability  tool  inputs  and  outputs. 


Case 

Pol 

Ship 

U 

[m/s] 

<p 

[degl 

MSVR 

Lks 

V 

Mini 

[ml 

Min  PCR 
TdBl 

1 

HH 

Le  Rif  Ferry 

1 

21 

1.72 

1.76 

119 

267.6 

14.1 

HV 

1.98 

1.76 

00 

85.4 

16.2 

2 

HH 

Le  Rif  Ferry 

4.9 

0 

1.94 

2.25 

20 

55.4 

15.0 

HV 

2.44 

2.25 

00 

54.9 

21.8 

3 

HH 

Le  Rif  Ferry 

6.3 

0 

3.75 

3.73 

00 

17.8 

13.2 

HV 

4.50 

3.73 

00 

39.4 

21.5 

4 

HH 

Quest ,  Research 

6.1 

0 

2.11 

2.25 

49 

52.2 

14.3 

VV 

2.28 

2.25 

00 

59.1 

13.5 

HH 

Grenfell,  CCG 

2.11 

2.25 

49 

52.2 

14.3 

VV 

2.28 

2.25 

00 

59.1 

13.5 

HH 

Sambro ,  CCG 

2.11 

2.25 

49 

52.2 

14.3 

VV 

2.28 

2.25 

00 

59.1 

13.5 

5 

VV 

Nakhchivan ,  Cargo 

5.6 

0 

1.66 

1.76 

45.8 

399.0 

14.8 

VH 

1.81 

1.76 

00 

112.7 

15.4 

VV 

Stav.  Prince ,  Tanker 

1.66 

1.76 

45.8 

399.0 

14.8 

VH 

1.81 

1.76 

00 

112.7 

15.4 

VV 

Saral,  Ferry 

1.66 

1.76 

45.8 

399.0 

14.8 

VH 

1.81 

1.76 

00 

112.7 

15.4 

6 

HH 

Leo  Glory ,  Tanker 

15.3 

90 

10.4 

3.50 

00 

50.8 

12.3 

HH 

Sandettie ,  Cargo 

15.3 

90 

10.4 

3.50 

00 

55.6 

12.3 

B.1  Case  1  -  ASAR  dual  polarization  HH+HV,  low  incidence 
angle 

These  cases  are  summarized  in  Figure  27.  In  this  figure,  as  well  as  the  next  three  figures,  the  top 
row  shows  the  overall  image  (left)  and  a  zoom  of  the  outlined  area  (right).  The  image  channels 
are  colour  coded  in  red,  green,  and  blue  (RGB)  as  indicated  in  the  caption.  The  next  rows  show 
zooms  on  the  ships  of  interest,  either  for  the  available  polarizations  or  for  ships  considered,  while 
the  next  rows  show  point  target  analyses.  The  bottom  row  shows  screen  shots  of  the  Ship 
Detectability  Tool.  The  constituent  images  in  the  case  studies  are  presented  here  in  small  format. 
They  are  available  along  with  other  supporting  information  at  full  resolution  in  the  final  report 
presentation.  The  ships  in  this  image  were  validated  using  coastal  AIS  data. 

In  Case  1  there  is  a  known  ship  at  an  incidence  angle  of  15.6°  that  is  not  visible  in  the  co¬ 
polarization  channel  but  is  visible  in  the  cross-polarization  channel.  The  minimum  detectable 
ship  lengths  are  268  m  and  85  m  for  HH  and  HV  respectively,  with  the  actual  ship  length  of  129 
m  being  between  the  two.  This  case  illustrates  improved  ship  detection  performance  at  smaller 
incidence  angles  for  cross-polarization  data. 
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B.2  Cases  -  ASAR  dual  polarization  HH+HV,  medium 
incidence  angle 

This  case  is  summarized  in  Figure  28  and  shows  the  same  ship  as  was  considered  in  Case  1  but  at 
an  incidence  angle  of  29°.  The  minimum  detectable  ship  lengths  are  about  55  m  for  both  HH  and 
HV.  The  observed  PCRs  are  17  dB  and  20  dB  for  HH  and  HV,  respectively.  This  case  illustrates 
ship  detection  performance  that  is  comparable  for  the  two  channels;  the  ship  of  interest  is  readily 
detectable  in  both  channels. 


B.3  Case  3  -  ASAR  dual  polarization  HH+HV,  high  incidence 
angle 

This  case  is  summarized  in  Figure  29  and  shows  the  same  ship  as  was  considered  in  Cases  1  and 
2,  but  at  an  incidence  angle  of  44°.  The  minimum  detectable  ship  lengths  are  about  18  m  and  40 
m  for  the  HH  and  HV  channels,  respectively.  The  observed  PCRs  are  3 1  dB  and  24  dB  for  HH 
and  HV,  respectively.  These  cases  illustrate  improved  ship  detection  performance  at  larger 
incidence  angles  for  co-polarization  data;  the  ship  of  interest  is  readily  detectable  in  both 
channels. 


B.4  Case  4  -  ASAR  dual  polarization  HH+W,  small  ships 

This  case  is  summarized  in  Figure  30  and  shows  three  known  ships  (74  m,  68  m,  and  16  m)  that 
were  imaged  during  a  controlled  ship  detection  trial  held  off  Halifax  in  2004.  The  minimum 
detectable  ship  lengths  are  52  m  and  59  m  for  HH  and  VV,  respectively.  The  two  larger  ships  are 
readily  visible  in  both  available  channels.  The  smaller  ship  is  marginally  visible  in  the  HV 
channel.  This  is  consistent  with  the  minimum  detectable  ship  length  relative  to  the  actual  vessel 
sizes. 

B.5  Case  5  -  ASAR  dual  polarization  W+VH 

This  case  is  summarized  in  figure  and  shows  three  known  ships  (244  m,  140  m,  and  119  m)  that 
were  validated  using  coastal  AIS  data.  The  minimum  detectable  ship  lengths  are  399  m  and  113 
m  for  VV  and  VH,  respectively.  With  a  PCR  of  15.8  dB,  even  the  largest  ship  is  barely  visible  in 
the  VV  channel.  All  three  ships  are  readily  visible  in  the  VH  channel,  which  is  consistent  with 
the  minimum  detectable  ship  length  relative  to  the  actual  vessel  sizes. 

B.6  Case  6  -  RADARSAT-1  ScanSAR  Narrow 

This  case  is  summarized  in  Error!  Reference  source  not  found,  and  shows  a  tanker  (334  m)  and 
a  cargo  ship  (142  m)  in  RADARSAT-1  SCNA  imagery.  The  minimum  detectable  ship  lengths 
are  51  m  and  56  m,  respectively,  and  both  ships  are  readily  detectable  with  PCRs  of  31  dB  and  25 
dB,  respectively.  These  cases  illustrate  good  ship  detection  performance  across  a  broad  image 
swath  of  300  km. 
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Figure  27:  Case  1  with  HH  (left)  and  HV  (right).  RGB  =  (HV,  HH,  HV).  (Continued  on  next 

page.) 
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3AR  Ship  Detectability  performance  TcoE 
Sensor-Spec  I  fie  Calculator 


[select  Beam  j 

Inoicrervc  e  ftnflle  [deg]  “"J : 
RKOa  [m] 

RiiOr  (GrndRngl  [m] 

H  E  SZ  [(10  J 
MLOOKS 


Wind  Speed  [m/s] 
li'Jind  Aspect  Angle  [deg] 
K'PI&trJtnJlEQH  Order  Parameter 
Probability  of  Falsa  Alarm 
Margin  [dB] 


P 


Mi nl m urn  Detectable  Ship  Lcngili  [m] 
Mlnlmyni  Detectable  rtes  [da-fnJ'2I 
Clutter  to  Noise  Roiio[d&l 
ocean  quitter  [da] 

Minimum  p*aK  To  Cluaer  Ratio  [dej 


Ctear 


Figure  27:  Concluded. 


Figure  28:  Case  2  with  HH  (left)  and  HV  (right).  RGB  =  (HV,  HHf  HV).  (Continued  on  next 

page.) 
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Wind  Speed  [mfs] 

mrm 

Mnd  Aspect  Angle  [deg.] 

K-aisETtUmton  Order  Parameter 

■  I'T'l'B 

Prob  ability  of  False  Alarm 

Margin  [dB] 

KB 

1  .Reset  DeTauELs 

Figure  28:  Concluded. 


Figure  29:  Case  3  with  HH  (left)  and  HV  (right).  RGB  =  (HV,  HH,  HV).  (Continued  on  next 

page.) 
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SAP  Ship  Detectability  performance  Toot 
Sensor-Spec  If  I  e  Calculator 


knvis*r-*sui  _ 

fsefee  t  Seam  | 


fnelBerKe  Annie  [Beg]  f 
SKOa  [m] 

RiLQr  (GrndRo^l  [mj 

HE'SZ[dBJ 

NJ.QCK3 


Wind  Spc  cd  [mfs] 

■ins 

Wind  Aspect  Angle  [Beg] 

KM 

K'Distrlbuiton  Order  parameter  ^LLifl 

Probability  cf  False  Alarm 

Margin  [dB] 

n 

i  ' 

—CaJcuIttcrd  Ro  suits— 


Minimum  Detectable  Ship  Lengm  |m]  | 
Minimum  Detectable  RCS  fd&kn^i] 
Ciuner  to  Noise  Retro  [dB] 
ocean  clutter  [di] 

Minimum  P*aK  TO  Cfufler  Ratio  [dBJ  ( 


f  Clear 


SAR  Ship  Detectability  performance  Tool 
Sensor-Specific  Calculator 


IHVlSfll'K 


[  sefect  Beam] , 


InciBens  e  Annie  [Beg]  ^ 

RKOa  Em] 

fiHQr  (GmdRng}  [m] 

N£SZ  fdBJ 
NLCCKS 


Wind  5pc  cd  [mts] 

Wind  Aspect  Angls  [Beg] 

K'Pi&trlBulEon  Order  parameter 
Prcb  ability  at  Falsa  Alarm  2E2E33 

Margin  [dB] 


Reset  Defaults 


-Cakutatfrd  Results^ 
Minimum  Ocicctabla  Ship  Lengila  |m]  | 
Minimum  Detectable  ttes  [ciei-mjsai 
duflerto  Noise  Retro  [d&I 
ocean  Clutter  [da] 

tWptmum  Pean  To  Clutter  Ratio  [BBJ  ] 


I  Clt,;,f  i 


Figure  29:  Concluded. 
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Figure  30:  Case  4  with  Grenfell  (left)  and  Quest  (right).  RGB 

next  page.) 


(VV,HH,  VV).  (Continued  on 
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SAR  Ship  Detectability  performance  Tool 
Sensor-Spec  Ifle  Calculator 


Serect  Bearn  | 


lncltrerKeftnflle[tteci] 
RKOa  [tn] 

RliOr  (GrndRngl  [m] 

M  E  SZ  [dB] 

BLOCKS 


Wind  Speed  [m/s] 
li'Jind  Aspect  Angle  [deg] 
K'GII&trJtnJlEnn  Order  Parameter 
Probability  of  Falsa  Alarm 
Margin  [dB] 


CdEcidatad  Rest 
Minimum  De Lt c [able  SNp  Length  |m] 
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Figure  30:  Concluded. 


Figure  31:  Case  5  with  VV  (left)  and  VH  (right).  RGB  =  (VV,  VH,  VV).  (Continued  on  next 

page.) 
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SAft  Ship  Detectability  Performance  Toot 
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Figure  31:  Concluded. 
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Figure  32:  Case  6  with  a  tanker  (left)  and  a  cargo  ship  (right).  (Continued  on  next  page.) 
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Figure  32:  Concluded. 
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List  of  acronyms 


AIS 

Automatic  Identification  System 

AMAP 

Analysis  of  Marine  Applications  Potential 

AP 

Alternating  Polarization  mode  (Envisat  ASAR) 

APP 

Alternating  Polarization  Product  (Envisat  ASAR) 

ASAR 

Advanced  SAR  (Envisat) 

CCG 

Canadian  Coast  Guard 

CNR 

Clutter-to-Noise  Ratio 

Co-Pol 

Co-polarization  (i.e.,  HH  or  VV) 

Cross-Pol 

Cross-Polarization  (i.e.,  HV  or  VH) 

DRDC 

Defence  Research  &  Development  Canada 

ECDF 

Empirical  Cumulative  Density  Function 

EH 

Extended  High  mode  (RADARS AT) 

EL 

Extended  Low  mode  (RADARSAT) 

ENL 

Equivalent  Number  of  Looks 

ESA 

European  Space  Agency 

EW 

Extra  Wide  swath  mode  (Sentinel- 1) 

F 

Fine  mode  (RADARSAT) 

FY 

First  Year 

GMES 

Global  Monitoring  for  Environment  and  Security 

HH 

Horizontal  transmit,  Horizontal  receive 

HV 

Horizontal  transmit,  Vertical  receive 

IMP 

Precision  Image  Product  (Envisat  ASAR) 

IS 

Image  Swath  mode  (Envisat  ASAR) 

IW 

Interferometric  Wide  swath  mode  (Sentinel- 1) 

MY 

Multi  Year 

NESZ 

Noise-Equivalent  Sigma-Zero 

PDF 

Probability  Density  Function 

PRF 

Pulse  Repetition  Frequency 

R1 

RADARSAT- 1 

R2 

RADARS  AT-2 
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R&D 

RCS 

RGB 

S 

s 

SAR 

SCN 

SCNA 

SCNB 

ss 

sew 

SCWA 

SCWB 

SNR 

VH 

VHF 

VV 

W 

W 

WSM 


Research  &  Development 
Radar  Cross  Section 
Red,  Green,  Blue 
Strip  map  mode  (Sentinel) 

Standard  mode  (RADARSAT) 
Synthetic  Aperture  Radar 
ScanSAR  Narrow  (RADARSAT) 
ScanSAR  Narrow  A  (RADARSAT) 
ScanSAR  Narrow  B  (RADARSAT) 
ScanSAR  Sub-Swath  (Envisat  ASAR) 
ScanSAR  Wide  (RADARSAT) 
ScanSAR  Wide  A  (RADARSAT) 
ScanSAR  Wide  B  (RADARSAT) 
Signal  to  Noise  Ratio 
Vertical  transmit,  Horizontal  receive 
Very  High  Frequency 
Vertical  transmit,  Vertical  receive 
Wide  mode  (RADARSAT) 

Wave  mode  (Sentinel- 1) 

Wide  Swath  Mode  (Envisat  ASAR) 
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